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region. They show that in the spin-1 ladder, a sufficiently strong uniform field engenders a c = 1 
critical state regardless of the sign of J±, whereas the staggered field is expected not to yield any 
singular phenomena. From the bosonization techniques, new field-theoretical expressions of string 
order parameters in the spin-1 systems are also proposed. 
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I. INTRODUCTION 



Spin ladder systems have been investigated theoreti- 
cally for more than a decade. Several real magnets corre- 
sponding to them have been synthesized and observed. 
A trigger of such studies may be the discovery of high- 
Tc materials and the connection between the Hubbard 
model (one of the high-Tc models) and the antiferromag- 
netic (AF) Heisenberg model. ^ Recently, it has been rec- 
ognized that spin ladders themselves can provide several 
theoretically interesting phenomena: for example, quan- 
tum critical phenomena, nontrivial magnetization pro- 
cesses (plateaux and cusps), the connection v^rith field 
theories and integrable models, topological or exotic or- 
ders, etc. Especially, the intensive studies have largely 
developed the physics of 2- leg spin-i ladders. 

In the spin ladder systems, like other magnetic sys- 
tems, their responses to external magnetic fields have re- 
ceived theoretical and experimental attentions. Recently, 
in addition to the standard uniform magnetic field, stag- 
gered fields, which have an alternating component along 
a direction of the system, have been in the spotlight. ^"^^ 
Actually, several phenomena induced by them have been 
observed, and some mechanisms generating them in real 
magnets are known. ^'^^ As will be discussed in the next 
section, uniform- or staggered-field effects in the 2-leg 
spin-i ladder with AF intrachain coupling have been un- 
derstood well. In the uniform-field magnetization pro- 
cess, a massless phase exists between the saturated state 
and a massive spin-liquid state. A staggered field yields 
a quantum phase transition. 

Here, the following natural question would arise: how 
the external fields influence high-spin or iV-leg ladder sys- 
tems? To answer this (partially), this study specifically 
addresses low-energy properties of 2-leg spin-1 ladders 
with external fields through the use of several analytical 



tools. Our main target is the simple Hamiltonian, 

H = J ^ Sij ■ Sij+i + J± ^ Sij ■ S2,j + Hz, (1) 

l,j 3 

where Sij is the spin-1 operator on the site {I, j); I (= 1 or 
2) is the chain-number index and the integer j runs along 
each chain. The intrachain coupling J is positive. We will 
refer to such ladders as "AF" ladders. The interchain 
coupling J± is called the rung coupling. The Zeeman 
term Hz here is chosen as two types 

Hu = -HJ2{SIj+SIj), (2a) 
Hs = -Hj2{-^y{Sl^+Sl^\ (2b) 

3 

where H{> 0) is the strength of external fields. The 
latter (2b) has a staggered field along the chain (— l)^i7. 
Obviously, the AF rung coupling and the staggered field 
compete with each other. 

Regarding the case without external fields, some pre- 
ceding theoretical studies^^"^* of the spin-1 ladder (1) 
exist (although, as will see in Sec. Ill, there are also a 
few studies dealing with Hu)- Their results deserve to 
be summarized here for our consideration in later sec- 
tions. These studies explain how the first-excitation gap 
varies in dependence upon the rung coupling J±. At the 
decoupled point J± = 0, the model (1) is reduced to 
two spin-1 AF chains. As known well, the chain has a 
finite first-excitation gap (Haldane gap).^^ Around the 
decoupled point, the gap reduction takes place with \J±\ 
increasing. Namely, the gap has a cusp structure around 
J± ~ 0. Far from the decoupled point, in the AF-rung 
side, the gap increases together with the growth of J±. 
It approaches the gap of the rung dimer (two spins along 
the rung) which is the strong AF rung-coupling limit of 
the model (1). On the other hand, for the ferromagnetic 



FIG. 1: Schematic gap profile in the spin-1 AF ladder. The 
symbol Anal denotes the Haldane gap (~ 0.41J). For details, 
see Ref. 18. 



(FM)-rung side, even away from the decoupled point, the 
gap decreases monotonically. It reaches the Haldane gap 
of the spin-2 AF chain with the bond J' = J/2, which 
is the strong FM rung-coupling limit. (In general, the 
2-leg spin-S" AF ladder with the intrachain coupling J is 
reduced to the spin-2S' AF chain with J' = J/2 in the 
strong FM rung-coupling limit. )^'''^^ The first-cxcitation- 
gap profile is summarized as Fig. 1. A recent quan- 
tum Monte Carlo analysis^* quantitatively surveys the 
ground-state (GS) properties of the model (1) without 
Zeeman terms. It estimates the first-excitation gap and 
the spin-spin correlation length. Moreover, it proposes a 
new string- type parameter, which is discussed in a later 
section. The authors conclude that the GS is always 
massive from Jj^ = to cx3, and is characterized by 
the "plaquette-singlet" solid (PSS) state, one of short- 
range resonating-valence-bond (RVB) states, which is 
explained in Fig. 2. It smoothly connects two limiting 
states: the Haldane state (Jj^ = 0) and the rung-dimer 
state (J_L — > oo). The new string parameter is able to 
capture a feature of the PSS state. 

In addition to the external- field effects, we will reex- 
amine and revisit some of the results without external 
fields. 

Our analysis specifically addresses two regions: the 
strong AF rung-coupling region ( Jj^ 3> J) and the weak 
rung-coupling one (J 3> |J_l|). For the former (latter) 
region, we mainly employ degenerate perturbation the- 
ory (DPT) (field-theoretical methods). In the former 
region, we can also investigate 2-leg spin-S" ladders. It 
contributes to a systematic understanding of spin ladder 
systems. Furthermore, the comparison of our results and 
existing results of spin-i ladders will help to elucidate 
the spin-ladder systems. 

The organization of the remainder of this paper is as 
follows. First, we give a brief review of 2-leg spin-^ lad- 
ders in Sec. II. As mentioned previously, it is useful 
in comparing our spin-1 case with the spin-i case and 
highlighting the features of the spin-1 case. Section III 
specifically addresses the strong AF rung-coupling region 
(J_L > J). For the spin-1 ladder (1), the DPT pro- 
vides effective Hamiltonians and predicts that there are 
two critical areas (points) in the magnetization process 
applying the uniform (staggered) field. We also apply 
the DPT to the 2-leg spin-S" ladders, to find that in the 
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FIG. 2: RVB pictures in the GS of the spin-1 AF ladder (1) 
without external fields. The black point denotes a spin-i 
state. The ellipse including two black points indicates the 
symmetrization of two spin-i states and recovers an original 
spin-1 site. The black line represents for the singlet bond, 
which consists of a spin-i pair. For < Jx < cxd, the GS is 
well described by the PSS state. With increasing the rung- 
coupling Jx, the weight a (6) increases (decreases) in the GS 
wave function, and it approaches the tensor product of singlet 
dimers [singlet of the spin-1 pair]. For details, see Ref. 18. 

spin-S" case there are 25" critical regions (points) through 
the uniform(staggered)-field magnetization process. In 
addition, we propose RVB-typc pictures of noncritical 
phases in the uniform-field case with Ti^. In Sec. IV, 
we consider the weak rung-coupling region (J ^ |J_l|) 
and employ some field-theoretical approaches. In par- 
ticular, to map ID spin-1 systems onto a field theory, 
we utilize a non-abelian bosonization (NAB),^"^"^^ i.e., a 
Wess-Zumino-Novikov-Witten (WZNW) model descrip- 
tion. The first two subsections are devoted to an ex- 
planation of the NAB. Through it, the spin-1 ladder (1) 
is described by using a fcrmion field theory, which was 
originally proposed by Tsvelik.^^ Subsequently, we will 
consider the uniform-field case in Sec. IV C. In this case, 
the NAB is considerably effective because it can treat the 
uniform Zeeman term nonperturbatively.^^ From conse- 
quences of the DPT in Sec. Ill, the NAB here, and the 
gap profile in Fig. 1, we can determine the whole GS 
phase diagram in Sec. IV D. The critical regime is "sim- 
ply connected" and has a c = 1 criticality except for the 
decoupled point Jj^ 0. In Sec. IV E, we devote our 
attention to spin-1 ladders without external fields. Sec- 
tion IV E 1 is assigned to the evaluation of string-type pa- 
rameters in one-dimensional (ID) spin-1 systems within 
our field-theoretical framework. We propose their field- 
theoretical expressions. Using a RG analysis, we con- 
sider the GS phase diagram of a spin-1 ladder extended 
from the model (1) in Sec. IV E 2. Field-theoretical ap- 
proaches used here are not sufficiently efficient for the 
staggered- field case with Tig. A brief discussion on it is 
in Sec. IV F. In Sec. V, we summarize all results and 
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briefly discuss them. The Appendix offers readers sup- 
plements of field-theoretical techniques and calculations 
in Sec. IV. 



II. REVIEW OF SPIN-1/2 LADDERS 

In this section, we review low-energy properties of the 
spin-i AF ladder which is equivalent to the model (1), 
in which spin-1 operators are replaced with spin-^ oper- 
ators. 

For the case without external fields, a finite rung cou- 
pling engenders a gapped spin-liquid (no long-range or- 
ders occur) GS irrespective of the sign of J±. This is 
true because the GS of the spin-i AF Heisenberg chain 
is critical (massless) and the rung coupling is relevant 
from the standpoint of the perturbative renormalization 
group (RG) picture. The spin liquid can be illustrated 
using the short-range RVB picture^*"^^ as in Fig. 3. The 
figure indicates that the excitation gap in the spin-liquid 
phase has the same order as the energy required to cut 
a singlet bond. In the AF-rung side, singlet bonds tend 
to occur along both chain and rung directions. In con- 
trast, two spins on the rung tend to construct a triplet 
state on the FM-rung side. The tendency removes singlet 
bonds along rungs, and allows those along the diagonal 
direction. The FM-rung spin liquid must connect with 
the GS of the spin-1 AF chain {J± — oo) smoothly. 
Therefore, we call it the Haldane phase. In Refs. 30 and 
32, the authors show that these two kinds of massive spin 
liquids [(a) and (b) in Fig. 3] can be detected by the two 
string-type parameters, 

= - lim (^;e-S:;i,V,5«^a\ (3b) 

\]-k\^oo \ I 

where a = x,y^ (' ' ' ) represents the expectation value 
of the GS. We defined two new operators 5"° = ^ + 
S^j and = S'lj + S'f j+i. Indeed, from Fig. 3, one 
can confirm that b^^^„ ^ and O^J^d = (C^dd ^ 
and Oovon = 0) are realized in the AF-rung spin-liquid 
(Haldane) phase. 

For the uniform-field case, low-energy properties have 
become well understood, but the quantitative GS phase 
diagram has not been constructed yet as far as we know. 
For details, see e.g., Refs. 33-37 and references therein. 
The schematic GS phase diagram is expected as in Fig. 4. 
The spin-liquid, Haldane and saturated phases corre- 
spond to massive plateau regions in the uniform-field 
magnetization process. The critical phase, except for 
the decoupled line J± — 0, can be regarded as an one- 
component Tomonaga-Luttinger liquid (TLL),'^^ which 
is identical to a conformal field theory (CFT)'^^ with 
the central charge c ~ I. In this area, the magneti- 
zation per rung (Sj) changes continuously. When the 
field is increased, quantum transitions take place at lower 



-1 +1 -1 +1 -1 +1 -1 

//////////// 

-1 +1 -1 +1 -1 

FIG. 3: RVB pictures of the GS of the spin-| ladder. Panels 
(a) and (b) are typical spin configurations in the AF-rung 
and FM-rung spin-liquid phases, respectively. The gray line 
denotes a singlet bond. Up and down arrows represent Sf j = 
+ \ and —\, respectively. Each number under the dashed 
loop encircling two sites shows the value of 5*1 = Si^ -1-5*2,^+1 
in (a), and one of Sj — Sfj + 51 j in (b). Removing all sites 
of = or Sj ~ 0, one can see a "hidden" Neel order 
(-1-1, —1, -1-1, —1, ■ ■ ■ )• For rtiore details, see Refs. 30 and 32. 



H 

2x (c=1 chain) 
2J 




saturated phase 
<Sj>=+1 

c=1 critical phase Hi 

'I'Si) = 
spin-liquid phase 





Haldane phase < Sj> = 



FIG. 4: Schematic GS phase diagram of the spin-i ladder 
with the uniform Zeeman term (2a). 



and upper critical fields {Hi and H2)- They are of 
a commensurate-incommensurate (C-IC) type.^''"*^ The 
upper critical field H2 can be determined by calculating 
the exact spin-wave excitation energy in the saturated 
(perfect ferromagnetic) state: H2 = 2J + J± for AF- 
rung side, and H2 — 2 J for the FM-rung side. The same 
logic shows that the upper critical field of the spin-1 AF 
chain with J' = J/2 (strong FM rung-coupling limit) is 
2J, which is consistent with the result H2 = 2 J. 

Wang et al. have investigated the spin-i staggered- 
field case with the term (2b) by employing the DPT 
and Abelian bosonization.'^*'^^''*^"^° They predict that 
the competition between the AF rung coupling and the 
field H creates a second-order quantum phase transition. 
It belongs to a Gaussian type with c = 1, and sepa- 
rates a Neel phase and the massive spin-liquid phase, 
which continuously connects with the spin-liquid (a) in 
Fig. 3. Furthermore, Ide, Nakamura, and Sato^^ recently 
determined the transition curve with high accuracy using 
the level-crossing method^^ and a new twisted operator 
method. ^'^ The curve starts from the origin in the space 
{J±,H) because both the rung-coupling term and the 
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FIG. 5: Schematic GS phase diagram of the spin-i lad- 
der with the staggered Zeeman term (2b). In weak and 
strong rung-coupling regions, the transition curve follows 
h ~ 0.60 X j]!^^ jj^l"^ and h J± - 0.50 x J, respectively, 
where the factors (0.60 and 0.50) are determined by our nu- 
merical calculations.^^ These results are consistent with ana- 
lytical ones in Ref. 14. 




H 



FIG. 6: Eigenstates and eigenenergies of Tlj with the strength 
of the uniform field H varying. The thick dotted line corre- 
sponds to the GS. Each vector j ■ ■ ■ , • • ■ ) represents a state 



staggered Zeeman term are relevant to the single chain. 
The GS phase diagram is given in Fig. 5. 



III. 



STRONG RUNG-COUPLING LIMIT 



This section employs the DPT for the strong rung- 
coupling region: Jj_ ^ J. Frequently in ID systems, 
theories of the strong coupling limit such as the DPT 
provide visualizations of GSs and low- lying excitations. 
This is the case in the strong rung-coupling limit of the 
model (1), as we discuss in this section. 



A. Spin-1 ladders 

First, we investigate the spin-1 AF-rung ladder with 
the uniform field. The DPT^^ for this model has already 
been performed, for example in Ref. 54. However, we 
reproduce it here as a preparation for later discussions. 

Under condition J± 3> J, the 0-th approximation takes 
the following Hamiltonian: 



^o-th ^ ^^^.^ 



(4a) 



n, = J^Si^j-S2,j-HiSl^+Sl^), (4b) 

where the original model is reduced to a set of two-body 
problems on each rung. Because the j-th rung Hamil- 
tonian T-lj has two conserved quantities, Sj and Sj, one 
can easily solve it. Each eigenstate has a one-to-one cor- 
respondence to a state |S,§^)j where S and §^ are mag- 
nitudes of the rung spin Sj and its z component, respec- 
tively. The solution is given in Fig. 6. The GS of the 
rung encounters two level crossings in the magnetization 
process. In parameter space {J±,H), we concentrate on 
vicinities of these two level-crossing lines; H = J± and 



H — 2J±. Near one line, H = 2J±, two low-lying states 
on the j-th rung are 

^|1,1)„ ^|2,2),. (5) 

Low-lying states near another line are 

|+)2,, ^|0,0)„ |-)2,, ^|1,1),. (6) 

The first-order calculation of the DPT is equivalent to 
projecting the total Hamiltonian (1) onto the subspace 
which consists of the set of two states (5) [or (6)] over 
all rungs. Two projection operators can be defined as 
A = Uj Kj and A - rij A,,, where 



Pi, 



- l+>2,,2j(+| + |->2j2j(-|. 

For convenience, we define new spin-i operators as 



(7) 



f^/ = l+)i,.i..(-l, u- = \-),,,,M- (8) 



Similarly, another spin-i operator Tj is defined by replac- 
ing the subscript (l,j) to (2, j) in Eq. (8). The relation 
between these pseudo-spin operators and original spin-1 
ones is 



P2,jSijP2,j 



+ | A,,5±.A,, = (-iyC/7,(9a) 



'2^ 



, P2,jSi jP2,j 



2(-l) 
V3 



■T^. (9b) 



We read that (Sj) = 0, 1 and 2 correspond to (TJ) — i, 
(T/) = -i (or ([//) = i) and (U^) = -i, respectively. 
Under these preparations, we can obtain the following 
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effective Hamiltonian near H = 2J±: 

^ufi = AwA = A,iWA,j 

j 

-H^^i ^ C// + const, (10) 

i 

where Ji = 2J, Ai = i and ifu,i = 2Jj_ - + | J. 
Similarly, near H = J±, we obtain 

j 

-H^,2^T? + const, (11) 

where J2 = § J, A2 = ^ and ffu,2 = J± - -ff + 4- Both 
(10) and (11) are a spin-i XXZ chain with a uniform 
field. ^'"^ Known exact results on the spin-| XXZ chain 
immediately lead to the following predictions. From the 
model (10), in the region |i?u,i| < ■Jii^ + Ai) the sys- 
tem has a c = 1 criticality, in which the uniform mag- 
netization and the critical exponents of pseudo-spin cor- 
relation functions vary continuously with i?u,i varying. 
Otherwise [i.e., |-ffu,i| > + Ai)], the magnetiza- 
tion is saturated. Similarly, from (11), another massless 
region with c = 1 exists in |iJu,2| < 1/2(1 + A2). By 
interpreting these in the context of the original ladder 
model in the uniform field, we predict that two c = 1 
critical areas exist: J± — < H < J± + and 
2J± - J < H < 2 J_L + 4 J. Therefore, in the strong AF 
rung-coupling region, the spin-1 AF ladder has an inter- 
mediate plateau region with (Sj) = 1 in addition to two 
trivial plateau regions: the saturated state and the PSS 
state described in Fig. 2. The presence of the interme- 
diate plateau is contrasted with the spin- 5 case (Fig. 4). 
Prom effective models (10) and (11), we also see that 
all critical phenomena in the magnetization process do 
not involve any spontaneous symmetry breakings. While 
effective models also predict the intermediate plateau 
vanishes at the point {J±,H) = {^J,^J), where the 
lower and upper boundary curves H = J± + and 
H = 2J± — J cross each other. However, this estima- 
tion is too rough because the lowest-order DPT is prob- 
ably valid only in the sufficiently strong rung-coupling 
cases. A recent numerical study^^'^^ evaluates the van- 
ishing point {J±,H) « (1.44 J, 2. 7J). It claims that in 
the subspace fixing the total magnetization, the transi- 
tion between the c = 1 region and the plateau is of a 
Beresinski-Kosterhz-Thouless (BKT) type.^^ 

Imitating the RVB picture in Fig. 2, we attempt to 
regard the intermediate plateau state as a state compris- 
ing bonds of nearest-neighbor spin-i pairs. One should 
note that the pairs contain triplet bonds as well as singlet 
ones because the plateau has (Sj) = 1. From Ref. 59, the 
necessary condition for a plateau state is 

Su — Mu = an integer, (12) 




<S3> <0> <0> <0> cO> 



FIG. 7: Intermediate plateau state in the spin-1 ladder (1) 
with the uniform field (2a): the case with J — » (a) and 
the case where J is finite (/?) . Thin black bonds mean singlet 
bonds of spin-i pair. The dotted bonds mean triplet bonds 

whore both z components of two spins arc +i. In the panel 
(/3), gray plaquettes represent a super-position state drawn in 
the lowest part of the figure. 

where Su and are, rcspcctivcdy. the; sum of spin mag- 
nitudes and magnetizations within the unit cell of the 
plateau state. According to this condition, relations (9), 
and two effective models, the plateau state should be in- 
variant under the one-site translation along the chain and 
the exchange of two chains. These allow us to illustrate 
the bond picture of the intermediate plateau as Fig. 7. 
The state (/3) connects with the state (a) smoothly by 
decreasing the weight b. Hereafter, we refer to bond pic- 
tures including singlet and triplet bonds as "RVB" pic- 
tures. 

Next, we turn to the staggered-ficld case, which has 
never been discussed in previous studies. In the 0-th 
order DPT, the even-j rung Hamiltonian has the same 
form as Eq. (4b) , but the odd-j one has the uniform field 
pointing to the opposite direction to Eq. (4b). Therefore, 
low-lying states in odd rungs must be modified as follows: 

|+)i,i=odd ^ |1,-1)„ |-)i,i=odd ^ |2,-2)„ (13a) 

|+)2,j=odd = |0,0)j, |-)2,j=odd ^ (13b) 

Consequently, projection and pseudo-spin operators are 
redefined: new pseudo-spins in odd-j sites satisfy the 
same- type relation as Eq. (9), in which —\Uj + | 
(-^T/ + i) and Uf {Tj) arc replaced with - \ 
(iT/ - \) and -Uf {-Tp), respectively. Using these 
new tools, we can also construct effective models for the 
staggcrcd-field cases in a similar way deriving (10) or 
(11). In the vicinity of the line H = 2J±, the effective 
Hamiltonian is 

where Up^ = (-1)J[/J'^ fj/ = [// and i?s,i = 2Jj_-H- 
|j. Similarly, near the line H = J±, the effective one 
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7i°^ is the same type as Til*^, in which (J7i, Ai,[/") 

(J2,A2,I7) and ^ H,^2 = J^- H-\j. Un- 
like the uniform-field cases, the above two models have 

effective staggered fields Hs.i and -ffs,2, respectively. Be- 
cause such alternating terms are relevant for the spin-i 
critical chain, infinitesimal values of staggered fields im- 
mediately yield a finite excitation gap. In other words, 
only when the field i?s,i(2) vanishes, the GS is critical. 
Therefore, there exist two critical lines (not areas) with 
c = 1: i7 = J_L - ij and 7? = 2J_l - fJ- The former 
(latter) critical hne satisfies (C//) = and (S'/) = |(-1)-' 
((f/) = and = ^(-1)^). Both lines are of a c = 1 
Gaussian- type transition in common with the spin-i case 
(Fig. 5). On the critical lines, the staggered susceptibil- 
ity {—^y-§H' diverges [see the next subsection]. Like 
the uniform- field case, any symmetry breakings do not 
occur at the transitions. Numerical^" and analytical^'®" 
works show that the staggered magnetization in 71!°^ (or 

2) changes continuously with the staggered field vary- 
ing. Thereby, we can confirm that the original staggered 
magnetization (— 1)-'(5|) has no plateau, contrary to the 
uniform-field cases. 

We should discuss effects of the higher-order terms of J 
in the DPT, although they have already been explained 
in Ref 14. Especially, let us consider whether any mech- 
anisms varying the properties of above c = 1 criticalities 
emerge or not, from higher-order effects. In the vicin- 
ity of each transition, both two low-lying states in the 
rung and the perturbative intrachain coupling part are 
invariant under spin rotations around the z axis of the 
total spin. Therefore, the U(l) symmetry can not be bro- 
ken by higher-order terms. For instance, an anisotropic 
XY exchange interaction, which brings a mass genera- 
tion, does not occur. The continuous U(l) symmetry is 
one of the characteristic natures in the c = 1 criticality. 
Because the original model (1) has a site-parity symme- 
try, we can also say that bond-alternating terms, which 
produce a mass too, do not appear. As long as we focus 
on the strong nmg-coupling cases (Jx ^ J), anisotropy 
parameter Ai(2) will not exceed the value of the BKT 
transition point Abkt = 1- These considerations imply 
that higher-order perturbations can not make c = 1 crit- 
icalities of transitions change, although they will modify 
parameters (Ji(2), Ai(2), i/u(s),i(2)) slightly and generate 
several new but small terms (e.g., next nearest neighbor 
interaction terms) in effective models. 

Using analytical tools, we can also evaluate sevciral 
physical quantities around each transition. We will dis- 
cuss them in the next subsection along with general spin- 
S cases. 



B. Generalization to spin-S ladders 

Viewing results of DPTs in spin-i cases (Figs. 4, 5, 
Refs. 14 and 34) and above spin-1 cases, we find it pos- 
sible to generalize them to 2-leg spin-^* AF ladders with 



external fields (2a) or (2b). Even for the spin-S' case, in 
which Sij means a spin-S operator, the rung Hamilto- 
nian is solvable through states |S, S^)j, in which S runs 
from to 25'. The number of independent states |§, S^)j 

is Em=o(2m + 1) = (2.5 + l)^ and the table of Clebsch- 
Gordan (CG) coefficients (or the Wigner-Eckart theo- 
rem) can determine the relation between them and the 
states of original two spins {Sij,S2j)- The eigen energy 
of the state |S, S^)j is 



J± 



§(§ -hi) - J±S{S +i)-m 



(15) 



where we consider the even-j rung. In a subspace with a 
fixed value of §, the GS of the rung is |§, although 
at = 0, all (2§ + 1) states are degenerate. Because 
the GS energy of the above subspace goes down with the 
slope — S when H is increased, the GS of the full space 
has 25 level crossings: H = J±,2Jj_, . . . , 25J_l. (This 
is a general result of Fig. 6.) Hence, one readily ex- 
pects that the spin-5 case has 25 critical phenomena in 
the uniform(staggered)-field magnetization process. The 
first-order DPT near each level-crossing line gives the ef- 
fective model, which is the same type as the spin-1 case: 
we always obtain a spin-| XXZ chain with a uniform 
(staggered) field for the uniform(staggered)-ficld cases. 
Therefore, even in the spin-5 cases, similar conclusions 
as the spin-1 cases can be available. In principle, one can 
obtain effective models in all vicinities of level-crossing 
lines. However, derivations of the models associated with 
lower-field crossing lines must calculate numerous CG co- 
efficients (syntheses of spins). Tables I and II depict only 
the effective models near the highest-field and second- 
highcst-ficld lines. Nevertheless, these two will be suf- 
ficient to conclude that the spin-5 case has 25 critical 
regions (points) in the uniform(staggered)-field magneti- 
zation process. Note that in these two tables, we use the 
same symbols (Ji(2), Ai(2)) -ffu(s),i(2)) in the spin-1 cases 
to represent efli'ective models. 

From these tables, we can extract the following infor- 
mation on the GS of the spin-5 case, (i) For uniform-field 
cases, the critical areas around the highest-field line and 
second-highest-field line, respectively, are 



25 Jj_ -J<H< 25 Jj_ + AS J, 
852 + 25-2 



(25- l)Jj_ - 
< (25 - 1) Jj_ 



45-1 

2452 -145-1-1 
45-1 



J <H 
J, 



(16a) 



(16b) 



which determine the width of the critical region F in the 
tables. The plateau regions exist outside the regions F, as 
with the spin-1 case, (ii) The larger the magnitude of spin 
5 becomes, the more the width F increases as a result of 
the growth of the effective coupling Ji{2)- The effective 
model approaches a spin-^ XY chain. Because plateaux 
are characteristic in "quantum" spin systems, the growth 
of F, or the decrease of plateau regions, means the system 
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TABLE I: Effective models around the highest-field lines H = 2SJ±. Parameters (»7i, Ai, J7u,i, J?s,i) are values of the effective 
exchange coupling constant, the anisotropy parameter, the effective uniform field and the effective staggered field, respectively. 
The sign T denotes the width of the c = 1 critical region along the H-axis in the original parameter space {J±,H), Mu,i 
represents the possible values of (Sfj) in the critical region F, and Ms,i is the value of KSf.j)! at the transition line Hs,i = 0. 



spin 


Ji 


Ai 






r 






Hb,1 




1/2 
1 


J 
2J 


1/2 
1/4 


J± 
2J± 


-H+U 
-H+^J 


3J 
5J 




J± ■ 
2Jx 


-H-U 
-H-U 
-H-U 

TT 4S-1 T 


1/4 
3/4 


3/2 


3J 


1/6 


3Jx 


-H + |,7 


7J 


1 — > 2 


3Jx 


5/4 


2 


4J 


1/8 


4Jx 


-H + lj 

TT . 4i'-l T 


9J 


2 — > 9 

2 ^ ^ 

o 1 . <:< 


4Jx 


7/4 
o 1 



S 25J l7(4S) 2S'Jx - + ^ J (i^Tlp S'- i 2S'Jx - -ff - ^ J 



TABLE II: Effective models around the second-highest-field lines Ji" = (25— l)Jx. Parameters (J72, A2, -H'u,2, T, Mu,2, -ffs,2, Ms,2) 
have the same roles as (J7i, Ai, i?u,i, F, Mu,i, i?s,i, Ms,i) in Table I, respectively. Note that the spin-i case does not possess 
the effective model. 

spin J2 A2 Hn,2 T Mu,2 Hs,2 M3,2 

~ Jj 3/16 J±-H+y 0~TT j^-H-U 174 

3/2 fj 5/48 2Jx-if + |J fJ 1^1 2Jx--H"-|J 3/4 

2 fJ 7/96 3Jx-ii + fJ if^J 1^1 3J±-H-U 5/4 

5/2 aoj 9/160 4Jx -H + lJ if-J | ^ 2 4J± - H - |j 7/4 

S '-^^J Tefe) (25-l)Jx-g + ^J ''':-s]_r-' J S-l^S-l (25 - 1) Jx - Jf - ^ J g-| 



approaches the classical vector spin system. On the other 
hand, in the case fixing S, the highest-field critical region 
r is smaller than the second highest-field one. Therefore 
we expect that the lower-field critical regions arc larger, 
(iii) Like the spin-1 case, all intermediate plateaux must 
vanish at a sufficiently weak rung-coupling region, (iv) 
For the staggcred-ficld cases, the transition lines with 
the highest field and the second-highest field are, respec- 
tively, 

H = 25Jx - (45 - 1) J/2, (17a) 
H = (2S'-l)Jx-(45-3)J/2. (17b) 

(v) All critical phenomena do not involve spontaneous 

symmetry breakings as with the spin-1 case. 

Now, let us investigate the physical quantities around 
each criticality in more detail. If we represent the effec- 
tive model near each transition with the (second) highest 
field using the pseudo-spin operator [7^" {T°') in a simi- 
lar manner of the preceding subsection, the original spins 
are projected out as follows: 



^, S^-^i-iy^/SUf, (18) 



where we consider even-j sites, again. For the staggered- 
field cases, as Eqs. (13) and (14), we must redefine J7j* 

(Tf ) in ah odd-j nuigs and define (Tf ). The relation 
(18), of course, contains Eq. (9), and provides the mag- 
netic relation between the effective and original models. 
As mentioned earlier, the low-energy properties of the 
effective models arc elucidated well by several analytical 
methods. Therefore, the combination of the relation (18) 



and such methods can provide accurate predictions of the 
original spin-S ladders. 

First we discuss the uniform-field cases. According to 

Bethe ansatz'''*'^-'^'^^ and Abclian bosonization,'^®''^^'^''"^'' 
the low-energy physics of the massless area |^fu,i(2)| ^ 
^1(2) (1 + ^1(2)) is governed by the free boson field theory 
with the spin-wave [massless excitation] velocity ?;i(2) and 
the compactification radius -Ri(2) ■ (We refer the reader to 
Appendix A for an explanation of Abelian bosonization.) 
Utilizing their knowledge, one can derive the uniform sus- 
ceptibility formula^'^^ 



Xl(2) 



dH^A 



or 



ao 



(27r)2i??(2)«i(2) 



, (19) 



where ao is the lattice constant. Because the effective 
field -ffu.i(2) varies linearly with the original field H (see 

Tables I and II), the relation Xi(2) = ^~§jf^ is real- 
ized for the magnetization process varying only H,^^ at 
least in the first-order DPT. One thus can regard all the 
behavior of Xi(2) as those of the original susceptibility, 
except for the difference of the factor i. The linear re- 
lation between the effective and original fields is often 
used below. The Bethe ansatz can determine the radius 
i?i(2) and the velocity 111(2) as functions of Ji(2), ^1(2) 
and i?u,i(2)- Especially for i^u,i(2) = 0, the radius and 
the velocity are represented analytically as 



Ri 



(2) 



^1(2) 



1 arccos Ai(2) 



1/2 



1 - A2 



(20) 



2 arccos A 1(2) 

Inserting Eq. (20) to Eq. (19), we see that when S goes 
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from one half (one) to oo, the pseudo-spin magnetization 
curve slope [i.e., the susceptibility (19)] at the midpoint of 
the (second) highest-field critical regime [i.e., at H^,i{2) = 
0] decreases monotonically from -y=^ — 0.184 x 

(0.0954 X J-i) to 2:;^ il^kj")- Fixing S, we also find 
that Xi(-ffu4 = 0) is larger than X2{Hu,2 = 0). For exam- 
ple, in the spin-1 case, Xi{Hu.i = 0) = 0.119 x and 
X2{Hu,2 = 0) = 0.0954 x J~^. These imply a reasonable 
fact that the larger F is, the smaller is the magnetization 
slope at -ffu.i(2) = 0. Subsequently, let us consider how 
magnetization approaches the value of the plateau (sat- 
uration). Without utilizing formula (19), studies of the 
C-IC transition^°"''^ have shown that near the saturation, 
the magnetization behaves as 

(C//) (or (T/)) ~ |J?- - ^^u,l(2)^/^ (21) 

where |^fu\(2)| is the critical value Ji(2){^ + ^1(2))-^^ 
Hence, in the magnetization process, the original magne- 
tization {Sij) also behaves, near each plateau, ~ \H'^'' — 

Hl^^"^, where H" is the critical field of each c = 1 re- 
gion. This power law is a universal property. That is, 
it is independent of the spin magnitude S and the level- 
crossing line we choose in the two-spin problem of the 
rung. The spin-wave analysis can exactly calculate the 
gap in the saturated state of each effective model. For 
l-f^u,i(2)l > l-^u'^i(2)l' t^*^ excitation gap is estimated as 
l-'^u, 1(2)1 — \^ui(2)\-^^ Again, translating this into the 
original model, we see that when H is moved just outside 
each c = 1 region, the gap grows as [if — H"\.^'^ 

Performing the same spin-wave analysis in the satu- 
rated state of the original spin-5 ladders, we can deter- 
mine the upper critical uniform field He- The field gives 
the boundary between the saturation with (Sj) = 2S and 
the c = 1 phase just under it. The result is 

Hc = 4:SJ + 2SJ±. (22) 

Surprisingly, perfectly agrees with the upper critical 
field derived from the effective model around the highest- 
field line [sec the region (16a)]. In other words. He is not 
modified by the higher-order perturbation effects of the 
DPT. For the FM-rung side, the spin-wave theory implies 
that 

He = AS J, (23) 

where He does not depend upon the rung coupling J±. 

Equations (22) and (23) arc a generalization of the criti- 
cal field H2 in the spin-i case (Fig. 4). 

Next, we shift our focus to the st agger ed-field cases. 
When the effective staggered field -ffs,i(2) is vanishing 
[i.e., the GS is massless], the bosonization translates the 
pseudo-spin operator to the following boson representa- 
tion: 



where (t>{x) [x = j x 0,0] and ^1(2) arc the boson field 
and a nonuniversal constant, respectively. Here we 
neglect the so-called Klein factor From Eq. (24), 
a finite i?s.i(2) leads to a perturbation term proportional 
to i?s,i(2) sin(0/i?x(2)) for the effective boson field theory, 
and it then becomes a sine- Gordon model. The vertex 
operator sin((/i/i?]^(2)) has the scaling dimension Xi(^2) = 
— and is always relevant (i.e., Xi(2) < 2). There- 

1(2) 

fore the scaling argument^ tells us that any small fields 
Hs,i{2) yield an excitation gap mi(2) and the pseudo-spin 
staggered magnetization as: 

mi(2) ~ |^fs,i(2)|'/''-"^(^)^ (25a) 
(t^/)(or(T/)) ~ (-l)J'|if,,i(2)ri(^)/[2-u^)l,(25b) 

where if i^s,i(2) < 0, we take the replacement (—1)-' 
(-1)^+1 in Eq. (25b). Substituting R^2) in Tables I and 
H for a:i(2), we find that the larger the value S becomes, 
the slower the growths of both the gap and the magneti- 
zation become. In other words, the singularity of the 

staggered susceptibility l^k^.i^o (or ^^\hs,2^o) 
decreases with S increasing. Particularly, in the limit 
5 ^ 00, mi(2) ~ |Fs,i(2)I and (or |(f/)|) ~ 

|-'^s,i(2)|, which means the staggered susceptibility does 
not diverge at the "transition" line ^^3,1(2) = 0. These 
can be again interpreted as a sign of the approach to 
the classical spin system. On the other hand, fixing S, 
one observes that when the fields i?s,i = i?s,2 are small 
enough, mi and are larger than m2 and |(TJ)|, 

respectively. We therefore anticipate that the transition 
with the higher field has a stronger singularity. In the 
spin-1 case, we have nii ~ |iJs,il°'^''*; "^2 ^ \Hs.2\^'^^^, 
|(C//)| ~ |ifs,ir-^^' and |(f/)| ~ |i?s,2r^°'. In common 
with the uniform-field cases, -ffs,i(2) has a linear relation 
with H. In order to translate all consequences into ones 
of the model (1) in the original staggered-field magneti- 
zation process, it is sufficient to replace -ffs,i(2) and (Uj) 
(or (f/)), respectively, with H and {S[.j) - (-l)^'Ms,i(2), 
where -Mg ^(2) is the staggered magnetization par site at 
each transition line (see Tables I and II). 

Summarizing all the above results about the spin-S* 
cases, we can draw GS phase diagrams and the uni- 
form and staggered magnetizations as in Fig. 8. Notably, 
Fig. 8 is valid in the strong AF rung-coupling limit. 

Utilizing solutions of the Bethe ansatz integral equa- 
tions,^ values of the nonuniversal constants ^1(2),^'*'^" 
etc., we can serve more quantitative predictions. We omit 
them here. 

Finally, we speculate the short-range "RVB" picture of 
the plateau states in the uniform-field case, without any 
computations. The spin-5 case has 25-1-1 plateau re- 
gions including two trivial plateaux: the saturated state 
and the H = state. The guides to guess the "RVB" 
pictures are the effective models, in which translational 
symmetry docs not break, the plateau condition (12), 
and the expectation that the bonds along the chain are 
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Uniform-Field Case 
(al) ^ saturation 

Ht 




(a2) 
H 



Staggered-Field Case 



c=1 curve 



2S 
2S-1 

3 
2 
1 




(pl) 



<s^,i-s:,> 



saturation 



2S-1/2 



2S-3/2 



3/2 



1/2 




(P2) 



■ transition point 

2SJi-(4S-1)J/2 
/ . 



(2S-1)Ji-(4S-3)J/2 



FIG. 8: Panels (al) and (a2) are, respectively, the GS phase 
diagrams in the uniform-field case (2a) and the staggered-field 
one (2b). Panels (/31) and (/32) are, respectively, the uniform 
magnetization curve in the case (2a) and the staggered one in 
the case (2b). In the panel (/32), we denote the magnetization 
of the even-j rung. 



subject to taking the singlet state and all plateaux vanish 
in the sufhciently weak AF rung-coupling region. From 
these, in order to build the "RVB" picture for the plateau 
with (5*1) = S in the spin-S case, we should perform just 

the following two procedures: (i) putting S triplet bonds 
and (25* — S) singlet bonds per a rung, and (ii) "joining" 
two nearest-neighbor rungs by singlet or triplet bonds, 
not to break the translational symmetry along the chain. 
The plaquette states of the lowest panels in Figs. 2 and 
7 are available for the procedure (ii). For instance, two 
plateau states with (S^j) = 1 and (5*1) = 2 in the spin-| 
case are described as (A) and (B) in Fig. 9, respectively. 
Similarly, the plateau with = could be captured 
by the set of a PSS state and a spin-hquid state (a) in 
Fig. 3. Following the above speculations, one would eas- 
ily produce "RVB" pictures for any plateaux. 

Within the "RVB" picture, every time the GS moves 
from a plateau to the plateau just above it with H in- 
creasing, one singlet bond along each rung is cut and 
exchanged for a triplet bond. This phenomena is rem- 
iniscent of successive transitions in the spin- 5* bond- 
alternated chain.^^''*^^'^^ 



IV. 



WEAK RUNG-COUPLING LIMIT 



This section describes the spin-1 ladder (1) under the 
opposite condition to the last section: J ^ |Jl|,7?. We 
start from the limit J± = H = 0. There, the ladder be- 
comes two decoupled spin-1 AF Heisenberg chains. The 
Haldane gap of the single chain complicates the map- 
ping to some field theories. Nevertheless, two famous 



(B) 



/i;,N 



FIG. 9: Expected "RVB" pictures for the plateaux with 
{S-) = 1 (A) and {S-) = 2 (B) in the spin-| case. 



mappings exist:^^ Haldane's method based on a non- 
linear sigma model (NLSM)^^ and the method applying 
a NAB.^^"^^ We will use the latter in this section. It is 
useful for treating several additional terms, in first prin- 
ciple, by the perturbation theory and the RG picture. 



A. Critical points of spin-1 AF Chain 

The application of the NAB to ID spin-1 sys- 
tems^^'^^'^'^'^'* is supported by low-energy properties of 
the spin-1 Takhtajan-Babujian (TB) chain^^'^^ and the 
following spin-1 bilinear-biquadratic chain: 



(26) 



where J > 0, Sj is the spin-1 operator on the site j, 
and the TB chain corresponds to 5 = —I ~ 6tb- This 
subsection present a brief review the model (26).^^ 

The efforts of several people^*"^^ since the late 1980's 
have advanced our understanding of the model (26). In 
particular, the low-energy properties in the region J > 
and 1^1 < 1 have been elucidated well. The GS phase 
diagram and the excitation gap in this region are sum- 
marized as in Fig. 10. At least three special points ex- 
ist aside from our target, the Heisenberg point 6 = 0. 
The TB chain is integrable and has massless excita- 
tions with the wave numbers (momenta) fc = and 
Tr/oQ. The point Jls = 1; called the Lai-Sutherland 
(LS) model, is also integrable. It has massless exci- 
tations with fc = and ±27r/3ao. The low-energy lim- 
its of TB and LS points are, respectively, equal to the 
level-2 SU(2) WZNW model^^-^^ cp^^^ 

1 SU(3) WZNW model. On the S axis, these two points 
are located in a quantum phase transition. The TB point 
separates the Haldane phase {\d\ < 1) and the massive 
dimerized phase {S < (5tb), which has twice degenerate 
GSs. On the other hand, the LS point separates the Hal- 
dane and massless "trimerized" phase {6 > Sls)- The 
transition belongs to a generalized SU(3) BKT type.^^ 
Features of the Haldane phase are the existence of a fi- 
nite gap between the unique GS and the isolated spin-1 
magnon mode, and a "hidden" AF order detected by a 
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massive 
dimerized 
phase 



gap 



massless 
trimerized 
phase 




Haldane phase 

FIG. 10: GS phase diagram and the schematic gap behavior 
in the model (26). 



string order parameter: 

0" = - lim /S*" e*'^^"=i+i^" S"?), 

^ {on, 



(27) 



which is the same form as the string parameter (3a) in the 
spin-i ladder, except that is a spin-1 operator. The 

AKLT^'* point (5akit = 5 is noteworthy because its GS is 
exactly identical to a valence-bond-solid state. Moreover, 
the point (5akit is the onset of an incommensurability: in 
(l >) 6 > (5akit, the real-space spin correlations has an in- 
commensurate fluctuation, that connects smoothly with 
the three-site period one at the LS point. 

Asides from points above, recent studies have de- 
scribed there are two characteristic points: (5dis — 0.4 
(Ref. 84) and Lifshitz point Suf = 0A4.^^'^^ In the re- 
gion — 1<S< Sdisi the momentum of the lowest magnon 
excitation stays at fc = n/ao. However, on the right side 
of (5dis, it has a deviation from k = tt/oo, and splits into 
two incommensurate momenta, which smoothly reach the 
massless points -|-37r/2ao and — 37r/2ao at the LS chain, 
respectively. At the Lifshitz point, an incommensurabil- 
ity appears in the spin structure factor. 



B. Effective field theory 

Reviewing the last paragraph, one notes that the TB 
point in the model (26) can be adopted as an underlying 
point to consider an effective field theory for the Heisen- 
berg chain {5 — 0). It is likely that the LS point is also 
available for the field theoretical description. However, 
as mentioned already, there are three points changing 
the low-energy properties between the LS and Heisen- 
berg points. Therefore, the LS point is not appropriate 
for deriving the effective theory of the Heisenberg point. 
This subsection provides the effective field theory for the 
spin-1 ladder (1) in our notation. 

The level-2 SU(2) WZNW model, which describes the 
TB point, has two primary fields: the 2x2 matrix field 
gmn [ni, n g {1,2}] with left and right conformal weights 
{j^, jg) and the 3x3 matrix field $mn [m,n G {1,2,3}] 
with weig hts (5,5). This WZNW model is identical to 



three copies of massless Majorana (real) fermions, as the 
WZNW model has c = 3/2 and the Majorana fermion 
theory, which is equivalent to a 2D critical Ising model, 
has c = 1/2. This identification allows a discussion of 
low-energy properties of the TB chain using the c = 1/2 
CFT, instead of the WZNW model. In imaginary time 
formalism where t = it {t is real time), we have the 
fermionic Euclidean action for the TB chain. 



Stb 



drdx vieLd.^L+CRd.CR], (28) 



where ^£(2;) and ^£(2:) [a — 1,2,3] are, respectively, the 
left mover of a Majorana fermion with weights (^,0) and 
the right mover with (0,^). The velocity of the exci- 
tations V is the order of Joq.^^'^^ Here, we introduced 
the "light-cone" coordinate: z — vt + ix, z — vt — ix, 

peated indices are summed. In the operator formalism, 
the fermions satisfy the following equal-time anticommu- 
tation relations: {^ai^),^^^!)} — ^ab^apSix — y). On the 
top of the action, there exist correspondences among the 
WZNW primary fields and fields of the c = 1/2 CFT. In 
accordance with Refs. 96, 74 and 17, $m„(z, z) comprises 
a bilinear form of ^£(z) and £,%{z), which indeed has the 
weights {h:h)- In addition, from Ref. 17, we have 



{z,z) (X ^{Ta)r. 
a=0 



iGa{z, z), 



(29) 



where (To)r 



and Ti 2 3 are the Pauli matrices. 



The fields Qa{z,z) are determined as 



(30) 



where aa{z,z) and iJLa(z,z) are, respectively, the order 
and disorder fields in the critical Ising model (c = 1/2 
CFT). Here, we use the cyclic index a + ?> = a. Both 
fields a a and /ia have weights [j^, j^)- An imaginary 
unit i is embedded in Eq. (30) to let the field g be a 
SU(2) matrix. The SU(2) current operators J1{z) and 
J'j^{z) in the WZNW model can be defined by fermions 
as follows: 

Jl{^^) = -^^aUUl: Jm^-^eaUUh (31) 

where eabc is the totally antisymmetric tensor and €123 — 
1. The currents JI^r{x) satisfy level-2 SU(2) Kac-Moody 
algebra. Through the NAB,^^'^^'^^ the spin operator in 
the TB chain is translated into the following sum of the 
uniform and staggered parts: 

-S^ « J£-HJS + iCo(-l)-'Tr[(.g-.gt)r,] 



ao 



Jl + JR + Cl{-lyaa^ia+l^ia+2, (32) 



where both Cq and Ci are nonuniversal constants. On 
the left-hand side, a = 1, 2, and 3 correspond to x, y 
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and z, respectively. This formula connects smoothly with 
the bosonized spin density of the spin-i ladder. From 
Eq. (32), one notes that the one-site translation causes 

Here we must mention a subtle point. The OPEs in 
Appendix B 2 show that a disorder field //o has an an- 
ticommuting character (so far we implicitly think of it 
as a bosonic object). One solution to maintain it and 
the Hermitian property of the staggered part of the spin 
density (32) is to modify the staggered part as 



(33) 



where the new parameter k has the same properties as 
an imaginary unit: n* = —k and = —1. We will 
sometimes use this modification below. 

Making full use of the above relations, one obtains low- 
energy properties of the TB chain. In order to achieve 
the field theory for the Heisenberg chain beyond the TB 
chain, one must add the following two perturbation terms 
to the action (28). 



(34) 



As long as the attached terms are restricted to rele- 
vant or marginal ones, only these two terms are admit- 
ted and possess spin rotational (see Appendix C) and 
one-site translational symmetries.®^ From the forms of 
Eqs. (32), (33), (34) and the Hamiltonian for the TB ac- 
tion (28), we can infer that the time reversal transforma- 
tion {S^,i) {-S^,-i) is mapped to {^l,^R,K,,i) 
{^[1,^1,— Similarly, we infer that link-parity 



transformation 



SI 



and site-parity one S°' 



Slj {Sg is fixed) correspond to [^l{x),^'^{x), g{x)] 
[TeRi-x),±eL{-x),9{-x)] and [a{x),Cni^),g{x)] ^ 
[TCfl(-a;), ±^2(-x), -g(-x)], respectively.®*^ 

Because the Heisenberg point {5 = 0) is far from the 
TB point (S = —1), m and A are phenomenological pa- 
rameters. It is known that one may take m > and 
A > in the Haldane phase (Fig. 10).^^ The inequahty 
A > means that the term — AJ£ is marginally irrel- 
evant. As shown in Eq. (B6), the Ising model picture 
tells us that m > indicates that each Ising model is in 
the disordered phase ^ 0. The mass parameter m 
must contribute to the Haldane gap. Consequently, three 
bands built of ^£ and can be regarded as the spin-1 
magnon modes in the chain (26). Running from S = —1 
to 5 = allows the velocity v to be rcnormalizcd. How- 
ever, it still has the order of Jao-'^^'^"" We will use the 
same symbol v for the renormalized velocity. In addition 
to V, for the other parameters (m. A, etc.), hereafter we 
use the same symbols whether they contain any renor- 
malization effects or not. It is widely believed that the 
formula (32) is applicable even for the Heisenberg chain 
because its low-energy excitations still stay around the 
uniform point k = and the staggered one n/ao (see the 
previous subsection). 

Heretofore, wc have obtained an effective field theory 
for the spin-1 Heisenberg chain. This framework was 



first proposed by Tsvclik^^ in 1990. Utihzing Eq. (32), 
we easily obtain the field theory for the spin-1 ladder (1) 
without external fields. The action is 

5'iad = Stb [Cl : ^r] + Stb [Cl 1 Cr] 

dTdx[zme^^% - XJIJ^ + imilCR - \JIJr 
+J^aoiW + + JlJl + JrJr) 

+Cfj±ao KKaalJ.a+ll^a+2^aij'a+lij'a+2 , (35) 

where quantities without (with) an overtilde ~ represent 

fields of the chain 1 (chain 2) with Z = 1 (Z = 2). We stress 
that the Hamiltonian for the action S'lad is invariant un- 
der spin rotation, one-site translation, time reversal, two 
parity transformations and exchanging the chain indices. 

On top of the rung-coupling term, it is possible to 
translate the other terms into field-theoretical expres- 
sions. The uniform Zceman term (2a) and the staggered 
one (2b) are, respectively, mapped to 

« -njdx {Jl + JI + JI + 4), (36a) 

w —CiH 1 dx {Ka3iiiH2 + K^3iJ'iP'2)- (36b) 



An advantage of the field theory used here is that the 
uniform Zeeman term is translated into the fermionic 
quadratic form (36a), which can be treated nonpertur- 
batively. Equation (36b) is not invariant under time re- 
versal and one-site translational operations. 

Other considerations regarding symmetries between 
the spin-1 ladder and its effective theory (35) are found 
in Ref. 17. 



C. Uniform- field case 

Sections IV A and IV B complete the main prepara- 
tion dealing with the spin-1 ladder (1). This paragraph 
presents a discussion of the ladder with a uniform Zee- 
man term (2a). Wc clarify what critical phase emerges 
when the uniform field is applied. It is easy to infer that 
a weak rung coupling does not collapse the Haldane gap 
of two decoupled spin-1 AF chains (see the Introduction 
and the next subsection). Moreover, in the single chain, a 
sufficiently strong uniform field engenders a c = 1 critical 
state. Therefore, we take the following strategy, 
(i) We review the effective theory for the c = 1 critical 
phase appearing in the spin-1 chain with a strong uni- 
form field. (ii) Then, adding the rung coupling terms 
perturbatively, we investigate the low-energy physics of 
the ladder (1) with a strong imiform field. 

Following the above scenario, first we explain how 
the c = 1 state is described within the field-theoretical 
scheme. Wc neglect the four-body interaction, the A 
term. One may interpret that it vanishes via the RG 
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procedure. Actually, it is believed that the effective the- 
ory without the A term is sufficient to describe the low- 
energy physics of the Heisenberg chain. In this case, 
the Hamiltonian for the chain 1 is 

Hchainl J dx ivilpldxi^L ~ i'lidxi^R) 

+*|(^i5.ei-a9.^^) + *meiel, (37) 

where, for convenience, we introduced a Dirac (complex) 
fermion. 

The Hamiltonian (37) is a quadratic form. The uniform 
field mixes two species ^ and ^£ ^j, while it docs not 
affect ^£ j^. Through Fourier transformations t^L.Rix) = 
/fc:an f e^'^V'L.fl(fc) anda,i^(^) = h^, ^{e^'''cL,R{k) + 
e"''"c[ ^(fc)} Ujix) = ei,fl,(a;)], and Bogoliubov 
transformations. 




where ZY(fc)ii,i2 = m/[2e(fc)(e(fc) T kv)Y'^, W(fc)2i.22 = 
±i[e{k) T H'/V[2e(fc)]'/', ^^3(^)11.12 = Z^(fc)ii.i2, 
Z^3(fc)2i,22 = i[e(fc)TM'/V[2e(fc)]'/' and e(fc) = [(fct;)2 + 
m^]^/^, we obtain the diagonalizcd Hamiltonian 

Hchainl = 1 ^ [e+{k)i^l{k)^+{k) 

+e_ (fc)vA (fc) V'- (fc) + e(fc)'^^ (fc)'^(fc)] + const, (40) 

where e±(fc) = e{k)±H. The three bands e(fc) and e±(fc) 
reproduce the Zeeman splitting of the spin-1 magnon 
modes. Their structures are given in Fig. 11. When the 
uniform field exceeds a critical value m, a Fermi surface 
appears in the lowest band e_ (fc) and its low-energy exci- 
tations then can be captured as a massless Dirac fermion. 
This corresponds precisely to the c = 1 state we are look- 
ing for. Here, the mass parameter m, after including all 
the renormalization effects, can be identified with the 
Haldane gap. 

Now, in this strong- field region {H > m), let us recover 
the A term, and take into account the rung coupling. 
Provided that we focus on low-energy physics, expanding 
the modes 4'+ (x) and V"- (x) , respectively, around fc = 
and the Fermi momentum k = kp = {H^ — m^)^/^/u will 




Dirac fermion 

FIG. 11: Band structures of the spin-1 AF chain with a uni- 
form field H. 

be allowed. Therefore, we obtain 

^Rix) « -Lv(a:) + e^'=--C/+i?t(-x) 

ipLix) « -^i!{x)-ie'^'="'U-R\-x) 

-ie-''""'U+L\~x), (41) 

where U± U{±kF)i2 = mly/2H[H + (ff2 ± rn?Y/% 
R(x) and L{x) are, respectively, the left and right 
movers of the Dirac fermion [%l)^{x) « e^^'''^ R{x) + 
Xji^x)]. In addition, we defined the field il){x) as 

/f e+(A:)vt(fc)^+(fc) « !dxf{^^^dl + m + i/)^. 
Similarly, we introduce fields L, R and -0 in chain 2. 
Substituting these new fields for the effective field the- 
ory (Tichaini + Hchain2 + A tcrm + rung couphng), we ob- 
tain the effective Hamiltonian (Dl) under the condition 
H > m. (Because it is too lengthy, the explicit form is 
put in Appendix D.) In the process, using Eqs. (41), we 
carefully approximated mass terms, interaction ones and 
currents J£ ^ while conserving the Hermitian property 
of each term. Moreover, we dumped all terms possessing 
some rapid fluctuation factors exp(±mfcFa:) {n is an inte- 
ger). The action toward the Hamiltonian (Dl) is written 
as 

SH>m = Si [L, R, L, R] + 52 , + 53 [V., 4'] 

+Sint[L, R, L, R] + S'fjjji, R, L, R; £,^,£.1] 
+Sl,[L,R,L,R;^,i:] 

+ SL , S ; V', ^] + 5fnt [V'^] + 5.tag , (42) 

where Sq, S^-^^ and S'stag denote the free part of each field, 
the four-body interaction terms and the staggered part of 
the rung-coupling term [i.e., the final term in Eq. (35)], 
respectively. Integrating out massive fields Co, "0) and 
-0 leads to the effective action containing only soft modes 
L, R, L, and R. Through a cumulant expansion, it can 
be expressed as 

Scff[L, R, L, R] ~ + S^^^ + {Sint)M 

where Sint = Sq=2 5'i7it + -S'stag, and (• • • )m indicates the 
expectation value of free parts of four massive fermions. 
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Because the Abelian bosonization is useful for interacting 
Dirac fermion models (see Appendix A), we introduce the 
boson fields 4> and (j) from the Dirac fermions [L, R) and 
(L, R), respectively. Applying the formula (BIO) and the 
known results of the 2-leg spin-i ladder with a uniform 
field [refer Eqs. (29) and (32) in Ref. 36] , we can bosonize 
the Ising-field product of 5stag as well as the fermion 
fields. As a result, Sstag is mapped as follows: 

Cfj±aoKk(Ta Ha+l Ma+25'a /io+1 /ia+2 ^ 

CfJ±ao{n3fi3 cos[y/T:{9 - 0)] + ^(Jsaa cos[^/^^{(|) - ^)]} 



+ fluctuating or irrelevant terms, 



where (0) is the dual of the boson field (p {(f), and the 
first (second) term on the right-hand side is generated 
from the x and y {z) components of the rung-coupling 
staggered part. The result (44) is also supported by the 
NLSM approach. ^"^'^"'^ The final fluctuating or irrelevant 
terms may be negligible in the low-energy limit. Prom 
Eqs. (42)-(44), up to the first cumulant, the effective ac- 
tion ScK yields the following bosonized Hamiltonian; 



H 



1st 



/ 



dx 



Anl + B+{dM^ + Cd^ct>+ (45) 

+Anl + B_{d^(j>_f 

+Ce_ cos(\/27r6'_) -|- cos{V^(t>-) , 



where we define the symmetric boson field = {(f) + 
4>) I \/2 and the antisymmetric one ^_ = (<?!>— ^)/ v^, and 
n± = dt(t>±/v' {9±) is the canonical conjugate (dual) of 
the boson (j)±. New parameters in Hist are 



A 

B± 

C 

Cfi,_ 



A 

47r 



v' A o o 
2 27r + " 

±^{ui+uir, 



{ul - u^f 



■y2 



Ce_ oc Ci Jj_ao(M3)M) 



(46b) 

(46c) 

(46d) 
(46e) 



where a in Eq. (46d) is the cut-off parameter in 

the bosonization formula (A3), 5^ = i{£,L£.n)M = 

■>'{^Ur)m > and 6^ = {i:ip^)M = {'>Pi'^)M > 0. 
Both (5| and 5^ are O(l/ao). Except for these two, 
the averages of products of massive fermions vanish in 
the first cumulant (>S'int)M- One should note {ij-3)m = 
{h)M 7^ and (c73)m = (^3)m = 0. In the deriva- 
tion of Eq. (45), for simplicity, we assumed that the 
Dirac fermion bands are always half-filled. From this. 



for example, we employed the relation L — LL^ = 
2L^L — 6o = 2 : L^L :, where the order of Sq is the 
inverse of the fermion wave-number cut off, and the sym- 
bol : : means the normal-ordered product. Observing the 
Hamiltonian (Dl) carefully and using the operator prod- 
uct expansion (OPE) in the c = 1 and c = 1/2 CFTs 
(Appendix B 2), we find that the second cumulant yields 
new interaction terms cos(-\/87r0_), cos(2-\/87r0_,L(K)) 
and cos[\/27r(<;^-,i:,(ij) + 3^!>_,H(t))] from ((S? J2)m and 
(S'?tS'stag)M- [Here, (j)±^L{R) is the left (right) mover 
of 4'± .] This means that the second cumulant does 
not produce any vertex operators having the symmet- 



(44) ric bosons 0+, 6+ and 



■LiR)- 



To verify this result, it 



is sufficient to note that the presence of exp(iv87r(/)+), 
exp(j\/87r^-l-), exp(i2\/87r(/>+,t) and exp(i2\/87r?!'+,_R) re- 
quires, respectively, fermion four-body terms L'^RL^R, 
L'^R'^Dr), L'^L'^DD and RRRR. It is expected that 
except for the above vertex terms, relevant or marginal 
interactions do not emerge in the higher cumulants. 

Besides the discussion related to the explicit count- 
ing of vertex operators, the Hamiltonian (Dl) and the 
cumulant expansion have the following four remarkable 
points, (i) Above vertex operators corresponding to some 
four-body fermions occur only through the rung coupling 
(meaning that the A term does not violate the c = 1 phase 
in the decoupled chain), (ii) We can apply an argument 
in Ref. 59, which cleverly employs the bosonization and 
symmetries of spin systems. From Appendix C, the U(l) 
transformation associated with the spin rotation around 
z axis is given by ipL,R — > tpL,Re-^'*' is a real num- 
ber) in the field theory, which accompanies 
i?t ^ i?te»v and V ^ ^e*^ via Eq. (41). Of course, 
fields with an overtilde ~ also receive the same transfor- 
mations. In the boson language, they correspond to a 
shift of the symmetric dual field 9^ ^ 9-^- + constant. 
This U(l) symmetry hence prohibits the emergence of 
any vertex operators with the dual field 9+. Similarly, 
let us consider the one-site translation. As one can see 
from Eq. (41), it causes L^{-x) L^~x - ao)e^*'=^'^« 
and R^—x) R^{—x — ao)e*'^*'"°. They are mapped to 
(f>+{x) (f)+{x + qq) + 2\/2kFao/r, where r is the com- 
pactification radius in the c = 1 theory considered now. 
It must be close to the value of the massless Dirac fermion 
l/\/47r (see Appendix A). Vertex operators with 0+ are 
therefore forbidden, except for those where kp becomes 
some special value. These prohibition rules are actually 
formed in the above counting of vertex operators. One 
can confirm that U(l) and one-site translational sym- 
metries are maintained in the Hamiltonian (Dl). (iii) 
Because the correlation functions of massive fields decay 
exponentially, no long-range interaction terms emerge in 
all the cumulants. For instance, the correlation lengths 
of two-point functions for or ip are at most the or- 
der of {mao)~^ X ao, where m has the order of the Hal- 
dane gap (= 0.41J),^°^ and mao (^ Joq) must be a con- 
stant in the present scaling limit. In fact, the correlation 
length of the spin-1 AF Heisenberg chain is only about six 
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times as large as gq.^'^^ The connected correlation func- 
tions of Ising fields 0-3 and fj,^ also decay exponentially. 
(iv) From (iii), roughly speaking, the expansion can be 
thought of as a Jl/ J expansion. 

From all the considerations below Eq. (45), when the 
rung coupling is sufficiently weak, the bosonized Hamil- 
tonian (45), possibly with vertex operators of antisym- 
metric fields only, could be adopted as an effective theory 
under a strong uniform field H{> m). Following a stan- 
dard prescription, we perform a Bogoliubov transforma- 
tion. 




(47) 



where K± — \/ A/B±, and the canonical relations are 
conserved: [4'±{x),U±{y)] = [(j)'±{x),U'^{y)] =i6{x-y). 
From the view of new boson fields (0j_ , W^. ) , interaction 
terms. 



cos(V8^0-) = cos{^/8^TK-(p'^), (48a) 
cos(V2^6i_) = cosU^9'_^, (48b) 



cos(\/87r0_) = cos 



(48c) 



have the scaling dimensions 2K-,1/{2K^) and 
respectively. The parameters J± and A yield a de- 
viation from K± = 1. It, however, would be quite 
small in the weak rung-cou pling re gion. The most rel- 
evant term, hence, is cos{ y/ 2Tr / K-6'_), and it locks the 
phase field d'_. Therefore, the (/)'_ mode always be- 
comes massive once the rung coupling Jl enters into 
the system. Other interactions cos(2V87r0__i(fl)) and 



/27r( 



,L{R) 



ML)) 



have a conformal spin, and 



cos[ 

it is reported that such fields may engender non-trivial 
effects. ^^'^"^^ However, they would not be as powerful as 
recovering any massless modes in the present case. 
On the other hand, for the symmetric part, the linear 
term Cdx4'+ is absorbed into the quadratic part by the 
shift 4>+ 4>+ + Cx/{2B+). From these results, we con- 
clude that for a strong uniform field {H > m), only the 
mode remains massless and a c = 1 phase is realized, 
irrespective of the sign of Ji . 

Finally, we note the limitations and the validity of the 
methods used here. The fermion band width in the effec- 
tive theory (37) is expected to be smaller than or equal to 
0{J). Hence, when the uniform field becomes too strong, 
such 'AS H > J, it is doubtful whether the theory (37) is 
valid or not. Furthermore, for such a strong-field case, it 
may be necessary to take into account other low-energy 
excitations. Meanwhile, the cumulant expansion is not 
reliable well when |Jj_| reaches 0(J). For the deriva- 
tion of Eq. (45) , we used the assumption that the Dirac 
fermion bands is half-filled. Removing it does not in- 
fluence the main results presented in this paragraph. It 
merely changes parameters in Eq. (45) a bit. 
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FIG. 12: Schematic GS phase diagram of the spin-1 AF ladder 
(1) with the uniform Zeeman term (2a). 



D. GS phase diagram of the uniform- field case 

The strategy of the preceding subsection is not suitable 
for determining the lower critical uniform field, while we 
have already known that the lowest-excitation-gap profile 
of the spin-1 ladder is given by Fig. 1. Furthermore, the 
lowest excitations must consist of a spin-1 magnon triplet 
[this expectation is trusted at least for the strong AF or 
FM rung-coupling regions. Moreover, from our effective 
theory (excitations of fermions ^£ ^ are interpreted as 
spin-1 magnon excitations) and the NLSM analysis, it 
would be also true for the weak rung-coupling region.] 
Therefore, the gap profile in Fig. 1 itself is equivalent to 
the shape of the lower critical uniform field in the space 
{J±,H). [As H is increased, one of the spin-1 magnon 
bands goes down linearly with H, a,s a. result of the Zee- 
man splitting. See the band e_(/c).] 

Taking into account the above lower critical field and 
predictions in Sees. Ill and IV C, we can finally draw 
the whole GS phase diagram of the spin-1 AF ladder (1) 
with the uniform field as in Fig. 12. The sharp form of 
the intermediate plateau area is one of the natures of 
the BKT transition: the correlation length outside of the 
critical phase in the BKT transition is still anomalously 
long, and it inversely indicates that the excitation gap, 
which is proportional to the width of the plateau in the 
present case, grows considerably slowly. Features of the 
spin-1 GS phase diagram are the phase boundary near 
the decoupled point and the existence of the intermediate 
plateau. On the other hand, the c = 1 universality in the 
critical phase is common to both spin-i and spin-1 cases. 



E. No-field case 

This subsection specifically addresses situations in 
which spin-1 ladders have no external fields. 
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1. String order parameters 

Here, we attempt to evaluate string-type order param- 
eters in spin-1 systems within the field-theoretical de- 
scription. 

First, we investigate that of the single spin-1 chain, 
Eq. (27). The estimation of the nonlocal part O^^ = 

exp(i7r Yln=]+i '^n) (note O^J = Og^) can be done simi- 
larly to the bosonization technique in Ref. 32, where the 

string parameters of spin-i ladders (3) were calculated. 

In the continuum limit, O^^ is approximated as 

« exp ( ± ZTT / dy[Jl{y) + {y)] ) , (49) 

J Xj +ao 

where x^jj = k{j) x ao and the staggered part of the spin 
density is dropped. Constructing the boson theory with 
the scalar field <p from two Ising systems {cra+ii(^a+2), we 
can translate Eq. (49) as 

« exp(T^^/^y dydyct>{y)) (50) 

~ e-xp[TiyfK{(j){xk) - <p{Xj))\ 

= {cOS.{^fK(j){Xk)) TiSin{y/TT(j){Xk))] 

X {cos{\fK(j){xj)) ± zsin(\/7r?!>(a;j))} 

~ {lla+\{Xk)iJ-a+2{Xk) T tJa+i (xfc)cra+2 (xfe)} 
X {lJ,a+i{Xj)Ha+2{Xj) ± aa+l{Xj)aa+2{Xj)} , 

where we used Eqs. (A4), (B8), and (BIO). The remain- 
ing problem is just two edge spins S*" and S^. [This 
problem does not appear in the calculation of Eqs. (3).] 
At least in the Haldane phase where {(Ja) = 0, the stag- 
gered parts of the edge spins could not contribute to O". 
While, the product of O^^ and edge-spin uniform parts 
can be evaluated using OPE rules: J£ ^ x aa+i(Ta+2 ~ 
and J£^^ x ^ia+iHa+2 ~ cra+i(Ta+2- Therefore, 
an expected field-theoretical form of is written as 
~ (Oe"x> ~ + (cTa+iCTa+a)^. The deriva- 

tion of this form, however, has some subtle aspects, which 
arc mainly attributable to the continuous-field (coarse- 
grained) scheme. A similar difficulty is also present in 
the estimation of Eqs. (3).^^'^^^ To eliminate it, some 
ideas that are independent of field theories arc necessary. 
Actually, Nakamura resolves such an ambiguity of the 
field-theoretical expressions of Eqs. (3) using a symme- 
try cleverly.^^^ For our spin-1 chain case, ^ in the 
Haldane phase, whereas = in the dimerized phase 
(sec Fig. 10). Counting on this fact, we can propose an 
appropriate form, 

0« - (Ma+lMa+2)', (51) 

which may imply that two edge spins follow the rule se- 
lecting the disorder-field portion from the exponential 
part Ogx- The formula (51) has the same form as the 
field-theoretical form of Eqs. (3). That similarity must 
be one reflection of the fact that the FM-rung spin-i AF 



ladder is tied with the spin-1 AF chain smoothly. Fur- 
thermore, it also reminds us that O" is exactly mapped 
into an FM order parameter by a nonlocal unitary trans- 
formation.^^^ 

Our new proposal (51) can also tell us the behavior 
of in the vicinity of the TB chain. A scaling argu- 
ment (or the exact solution for the 2D Ising model) leads 
to m (i5 -|- 1) and (na) ~ (^ + 1)^/^ near the chain. 
Therefore, we predict that the critical behavior, 

O" ^{S + l)^/^ (52) 

occurs in the Haldane phase close to the TB chain. 

We next examine the spin-1 AF ladder. We denote two 
string order parameters of chains 1 and 2 as and O2, 
respectively. In Ref. 18, the quantum Monte Carlo sim- 
ulation shows that: (i) a new string parameter (£'"£'2) 
is always finite for the AF-rung side; (ii) the string order 
parameter of each single chain vanishes, once an AF rung 
coupling is attached in the system; and (iii) in the AF- 
rung side, (0^02) decreases mitil J± ^ OAJ and then 
grows monotonically until J± 00 (the rung dimer) 
(see Fig. 6 in Ref. 18). From these results, it is expected 
that the new string parameter is a quantity character- 
izing the PSS state. We discuss how the field theories 
reproduce these, and what they can predict. Supposing 
that Eq. (51) is applicable even for the spin-1 ladder, we 
have 

{6^6^) ~ {Ha+lfl'a+llJ-a+2fj-a+2f 

(C0s(v^(/)a+l)c0s(v^(^a-h2))^ (53a) 

~ (cos(V^$)cos(V^$))2. (53b) 

In Eq. (53a), the boson (f>a is made from the two Ising 
systems (cTo, da)- In Eq. (53b), bosons $ and $ are made 
from (cra+i,(Ta+2) and (ffa+i, 0-0+2), respectively. Equa- 
tions (53a) and (53b) can be evaluated by a bosonized 
effective theory (El) [or (E2)] and another theory (E3) 
plus (E4); respectively. The semiclassical analysis for 
Eqs. (El) or (E2) predicts that <^a is locked to the point 
(j)a = for the single-chain case (Jj_ = 0), and at that 
time Eq. (53a) has a finite value. This is in agreement 
with the fact that {O^O^) = {Ol)"^ ^ is realized 
at the decoupled point. While, for a weak (but finite) 
rung-coupling case, the effective theory (E2) possesses a 
new potential proportional to J±sm{^(f)a), as well as 
the mass potential proportional to mcos(v^47r0a). Their 
combination must vary the locking point from <i>a = 
to a finite and small value irrespective of the sign of J± . 
Therefore, we reach the same conclusion as the first con- 
tent of (iii), and predict that the decrease of (OJOf) 
also occurs in the FM-rung side: the new string param- 
eter would have a cusp structure like the gap in Fig. 1. 
Similarly, let us also perform the semiclassical analysis 
for another theory, (E3) plus (E4). It predicts that $, 
and (pa are all locked to zero at the decoupled point. 
The rung coupling engenders a new potential (E4) pro- 
portional to J_LC0s[7r(0 — 6)], in which 9 and 6 are 
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dual fields of <I> and <I>, respectively. This potential tends 
to fix Q and G) instead of $ and $. Fixing 8 means 
that the fluctuation of $ is large because $ and 8 are 
a canonical conjugate pair (see Appendix A). Therefore, 
assuming that 8 — 8 is locked in the low-energy limit as 
the rung coupling is finite, we can predict that the rung 
coupling makes the string parameter of the single chain 
(d)") ~ (cos(-y7r$))^ become zero. This result consistent 
with the content (ii). Moreover, it implies that {O") also 
vanishes in the FM-rung side. Note that fixing 8 — 8 
does not mean (OfOf ) because it can be rewritten 
as (OfOf ) ~ (cos[\/^($ + $)]) + (cos[^/^($ - !>)]) and 
(cos[y^($ + $)]) is expected not to have a large fluctu- 
ation. 

We see that spin-1 and spin-| string parameters take 

quite similar field-theoretical expressions. However, one 
should bear in mind that the effective theory for the spin- 
^ ladder is different from that of the spin-1 ladder: the 
former is two coupled sine-Gordon-like chains;^^ the lat- 
ter is three coupled sine-Gordon- like chains in Eq. (El). 



2. In the vicinity of the TB point 

In this paragraph, we briefly consider the spin-1 lad- 
der, in which two spin-1 chains are located near the TB 
point (Fig. 10), through the perturbative RG technique. 
In the effective theory of such a chain, parameters m and 
A are much smaller than those of the spin-1 Heisenberg 
chain. Furthermore, as mentioned already, m ~ (5-1-1) 
is realized. For the region where m. A, and J± are con- 
siderably smaller than J, the perturbative RG method 
based on the TB fixed point becomes a reliable tool to 
investigate the low-energy physics. 

We construct one-loop RG equations for coupling con- 



stants applying the OPE technique. 
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First, we 



consider all relevant and marginal terms around the fixed 

point (28). They are summarized in Table III, where we 
classify them into nine operators {Oj : j — 1, . . . , 9} to 
render each operator invariant under the spin rotational 
transformation (see Appendix D) and the interchange of 
two chains. Moreover, we introduced energy operator 
(mass term) ea{z, z) — Operators 02.3,4.9 are gen- 

erated dynamically in the RG flow, even though they are 
not present initially in the action (35). 

In the low-energy effective action, the diinensionlcss 
coupling constants {Gj} toward the operators {Oj} can 
be defined as 



9 r 
Slad = S* + ^2 f 



dx vdr 



Tran 



(54) 



where 5* = Stb[^1,^%] + STB[^l,eR] is the fixed-point 

action, ckq O(ao)] is the short-distance cut-off param- 
eter, and Xj is the scaling dimension of Oj. The RG 



equations for {Gj} are the following: 

Gi = Gi-^(G2 + 4G3-2G5)Gi 

+tt{GI-GI), (55a) 

G2 = -2G? + Ag2G5 + ^(2G5-G7)G3 

+ 4G4G7 - Stt^GI - n^Gl (55b) 
^(3G7 - 2G5)G3 - ^G4G7 + 27r2G^, (55c) 



G4 



(2G2 + G3)G7 



1 

2^ 

+ 4(2G5 + G7)G4 - 27r2G8G9, 



(55d) 



G5 
Gr 



2G\ H — 12G2G3 — n^iGg — Gg) 



+ - 



.2 \ 2^2 



Gl + Gl + Gl 



:Gt 



^ (2G2 + G4) Gs + ^(3G^ + Gl + Gj) 

+27T^ (G8 + G9)G8, 



Ga — -Gs 



1 ^ ^ , G7G9 

^G4G9 + ^ 



(55e) 
(55f) 

(55g) 
(55h) 



ttGi 



1 



:G2 



^Ga - iG4 



1 



G7 Gs, 



Go = 



1 



SttGi + ^G2 + ^Ga - ^Gs ) G9. (55i) 



Therein, Gj = dGj (L) / dL and L is the scaling parameter 
(the infinitesimal scaling transformation is ao — * aoe'^^). 
We adopted a simple circle- type cut off.^^° 

Using the RG equations, we discuss the low-energy 
properties near two decoupled TB chains. Close to this 
fixed point sufficiently, we can approximate them as 



Gl 
Gs 
Go 



Gl 



n(Gl 



Gl), 



^ + -Gi)G8 

4 TT 



G9. 



(56a) 
(56b) 

(56c) 



Couplings Gl, Gs and Gg are more relevant than the 
other couplings. The coupling Gi bears a single-chain 
character, whereas both Gs and Gg are the representa- 
tives of the rung coupling. Moreover, we may omit Gg 
because its initial value is zero. Under these approxi- 
mations, the treatment of Eqs. (56) becomes fairly eas- 
ier. One can find two nontrivial fixed points (Gi,G8) = 
(— 57r/4, ±-\/5/4). Let us assume the presence of these 
two fixed points, even though they are not close to the 
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TABLE III: Operators in the RG procedure. The sign G-j[L = 0) means the initial value of each running coupling constant in 
the RG flow. The parameter ao is the cut-off parameter defined in Eq. (54). 

Operators Scaling dimension Xj Gj{L — 0) 

02 = (ELi£«)(ELiei>) 2 

03 = E!>6 2(eaei, + e.£a) + (e£C« + eieS)(f£lfl + lilS) 2 

04 = rn'R + CRdmiiR + inil) 2 

05 = JrJI + JrJI 2 -A X ^ 
Oe = JtJl + J-hJR 2 Jxao x f 
O7 = J^Jl + JITr 2 Jxao x^^^ 

On = Kkaa^la + lfJ.a + 2^afia + lfj.a + 2 3/4 Ci JxAO X Ofg X ^ 

Og = «:K(TiCr2Cr3(Ti(T2Cra 3/4 



TB point (Gi, Gg) = (0, 0). Linearization of the approx- 
imated RG equations around new fixed points indicates 
that both points are of a divergent type, thereby imply- 
ing the existence of phase transitions. Similarly, the TB 
point is of a divergent type. Therefore, there exist two 
phase transition curves connecting the TB point and a 
new fixed point or another. In the vicinity of the TB 



point where approximations Gi ~ Gi and Gs 



rGs 



are allowed, a conservation law under the RG transfor- 



mation, |Gi/Gj 



4/5, 



a constant, is realized. Taking into 



account this law and recalling that Gi oc m, Gg oc J_l, 
and TO ^ (<5+l) are realized near the TB point, we expect 
phase transition curves follow: 



(57) 



around the TB point. Consequently, we can draw the GS 
phase diagram near the two TB chains as in Fig. 13. In 
this figure, the right side S > —1 of the horizontal (decou- 
pled) line Jl = is probably not corresponding to any 
phase transitions. In fact, we know that on the Heisen- 
berg line 6 = 0, the point J± = does not correspond to 
any transitions. Although on the Heisenberg line, the GS 
of the strong AF-rung limit (rung dimer) is quite different 
from one of the FM-rung limit (spin-2 AF chain), both 
AF- and FM-rung sides may belong to the same phase. 
Whereas, it is not sure whether the left side 5 < —1 of the 
line Jl = corresponds to a phase transition or not. The 
dimerized phases 1 and 2 must break the one-site transla- 
tional symmetry along the chain direction. According to 
the Zamolodchikov's "c theorem," two critical curves 
starting from the TB point (c = 2 x 3/2), belong to a 
universality class with c < 3. 

From Fig. 13, it is believed that the area characterized 
by the PSS-state picture (or connected to a spin-2 AF 
chain) widely expands around two decoupled Heisenberg 
chains, in the space (6, J_l). 



F. Staggered-fleld case 

The low-energy action for the staggered field case is 
jiven by Eq. (35) plus Eq. (36b). The latter term is only 



dimerized 
plnase 1 



dimerized 
pliase 2 




(PSS state) 

8=-l 



TB point 

c=2x3/2 



5=0 



He isenberg 
ooint 



(spin-2 AF chain) 



FIG. 13: Expected GS phase diagram around two decoupled 
TB chains (c = 3). 



invariant under the U(l) rotation around the spin z axis; 
and it does not possess the SU(2) symmetry. This partial 
violation of the symmetry makes each operator Oj in Ta- 
ble III decoupled to two U(l) invariant parts through the 
OPE between Oj and the staggered term (36b). Neither 
part is invariant under the SU(2) rotation. We have to 
consider more than 21 coupling constants, to construct 
the RG equations. Although we actually constructed the 
RG, we do not record them here. We were unable to 
extract characteristic contents from them because they 
are extremely complicated. Coupling constants for mass, 
rung-coupling and staggered Zeeman terms all grow up 
just monotonically. Therefore, the GS in the staggered- 
field case must have some massive excitations. More- 
over, two critical curves in the AF-rung side (Fig. 8) do 
not reach the origin {J±,H) = (0,0). We infer that the 
present field-theoretical strategy can not predict how the 
curves finish. 



V. SUMMARY AND DISCUSSIONS 

We explored the spin-1 AF ladder (1) and some of its 
extensions. In the strong AF rung-coupling region, we 
performed the DPT and found that the GS phase dia- 
gram of the uniform(staggered)-field case includes two 
c = 1 critical areas (curves). Subsequently, we extended 
these results to the spin-S' ladders, and predicted that the 
spin-5 uniform(staggered)-field case has 25' c = 1 areas 
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(curves). Figure 8 summarizes GS phase diagrams and 
the magnetization curves. The upper critical uniform 
fields (22) and (23) were determined by the spin-wave 
analysis; we saw that, surprisingly, Eq. (22) is equal to 
predictions of first-order DPTs (see Table I). We also 
proposed the "RVB" picture of each plateau (massive) 
state. 



(S = an intege0^^^^3^^^3^^ [S = a half integer) ^^^.^J^^si 

H f : 



saturated 
phase 

4SJ 




saturated 
phase 

4SJ 




On the other hand, we applied field-theoretical meth- 
ods for the weak rung-coupling region. For the uniform- 
field case, we employed the non-abelian bosonization effi- 
ciently. Combining the consequences of weak and strong 
rung-coupling analyses, we complete the GS phase dia- 
gram of the uniform-field case as in Fig. 12. Meanwhile, 
our field-theoretical approach was not efficient for the 
staggered- field case. Therefore, the GS phase diagram is 
not determined perfectly. The uncertain area, i.e., the 
area where the critical curves in Fig. 8 (2a) vanish, is 
probably located in an intermediate AF rung-coupling 
region, which must be distant from both weak and strong 
rung-coupling regions. 

Using those field theories, we revisited and discovered 
some properties of ID spin-1 systems without external 
fields. We showed that field theories can describe string 
order parameters in spin-1 systems, and proposed for- 
mulas (51) and (53). We hope that these formulas are 
useful in the search for some new string-type parame- 
ters in spin-1 systems. We also considered the GS phase 
diagram around two decoupled TB chains (Fig. 13). 

Through the present work, one obtains GS phase di- 
agrams of both the spin-i and spin-1 ladders with the 
uniform Zeeman term (2a) (see Figs. 4 and 12). The 
former (latter) model consists of two gapless (gapped) 
spin chains. These and our prediction of Fig. 8 enable us 
to expect that phase diagrams of 2-leg integer-spin and 
half-integer-spin AF ladders are written as in Fig. 14. On 
the other hand, for the staggered-field case, we can give 
only the following two predictions about the spin-S" lad- 
ders. For the 2-leg integer-spin ladders, critical curves in 
Fig. 8 {2a) vanish in a weak AF rung-coupling region. 
For the half-integer-spin ladders. Figs. 5 and 8 imply 
that the origin {J±,H) = (0,0) is perhaps a multicrit- 
ical point, from which 25 critical curves start. 

Except for the spin-i cases, the bosonization tech- 
niques for higher spin systems have some subtle and phe- 
nomenological aspects. Establishing more sophisticated 
bosonizations is an interesting but difficult problem that 
remains for a future work. 

Nowadays, several spin-i ladder compounds have been 
reported. Unfortunately, the materials regarded as a 
spin-1 ladder have never been found. An organic com- 
pound "BIP-TENO"^ may be a candidate of spin-1 lad- 
ders, but its magnetic behavior differs from the simple 
ladder (1). Suitable theoretical predictions regarding it 
have not been constructed. 



FIG. 14: Expected GS phase diagrams of 2-leg spin-5' AF 
ladder with the uniform Zeeman term (2a). 
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APPENDIX A: ABELIAN BOSONIZATION RULE 

In this appendix, we briefly summarize the Abelian 
bosonization^^''^^''*'*"^° used in Sees. Ill and IV, which is 
deeply related with the concept of TLL and c = 1 GET. 
(Our notation is similar to Rcfs. 38 and 48.) 

Bosonization shows that (1+1)D Dirac fcrmion mod- 
els are equivalent to a (H-l)D boson field theory. The 
main root of this technique lies in the identification be- 
tween the massless Dirac fermion and the massless free 
scalar boson field theories. The former Hamiltonian is 
represented as 



dxiv{'ip\dxipL 



(Al) 



where and tpR^z) are, respectively, the left and 

right moving components of the Dirac fermion. The sign 
V denotes the "light" velocity. [In real time formalism, z 
and z means, respectively, i{vt + x) and i{vt — x).] The 
fermions obey the equal-time anticommutation relations 
{tpU^),ipl3{y)} = SapSix - y) and (x), i/'/3(2/)} = 0. 
The corresponding massless boson theory has the follow- 
ing Hamiltonian (here we do not make the terms of the 
zero-mode excitations^'*'*^ clear): 



cala 



dx - 

2 L 



(9,0)^ 



(A2) 



where (/)(z,z) — (j)L{z) -I- (j)R{z) is the scalar boson field, 
4>L(R) is the left (right) mover of (j) f^^id H = dt4>/v is 
the canonical conjugate of 0. The Hamiltonian (A2) 
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can be mapped to the same form for the dual field 
0{z,z) = 4>L ~ 4'R- The equal-time commutation rela- 
tions among boson fields are defined as [0ij(x), (j)R{y)\ = 
-[</>!, (a;),^L(y)] = |sgn(x - and [(t)L{x),(t)R{y)] = 0. 
Since above two theories have a chiral U(l) symmetry, 
this method is called the "abelian" bosonization. 
They have the following operator identities: 



Vl 



V2ttc 
VR 
V2 



: exp 



: exp 



-iV ^TTcjiLiz) 
iV47T(j)R{z) 



(A3a) 
(A3b) 



where r]L,R, called Klein factors, are necessary to 
guarantee the anticommutation relation between the left 
and right movers of the fermions, and they hence satisfy 
{Va, Vb} = "^Sab- The factor a in Eqs. (A3) is a parameter 
which order is the inverse of the wave-number cut off of 
the Dirac fermion. The exponential-type operators are 
called "vertex operators." Using a point-splitting tech- 
nique,^^''*® one can also bosonize the U(l) currents Jl 
and Jr as follows: 



Jl{z) = : ViV'L : {z) = -^d^(t)L, 

Jr{z) = : VjjV'i? : (z) = ^a,<^fl, (A4b) 

where the symbol : : stands for the normal-ordered prod- 
uct. [Note that the above U(l) currents are different 
from the SU(2) currents J£ and in Sec. IV. How- 
ever, if we think of the c = 1 CFT here as the c = 1 
part of the c = 3/2 WZNW model, one component of 
the SU(2) currents is proportional to the U(l) current. 
See Eq. (B8).] 

Prom these relations, one sees that Dirac fermion mod- 
els even involving arbitrary interactions can be mapped 
to a boson theory with vertex operators. The RG flow 
lets several kinds of interacting Dirac fermion models go 
to a free boson theory with a modified velocity [it is dif- 
ferent from V in Hamiltonians (Al) and (A2)] and a com- 
pactification radius R [which determines the period of the 
boson as (j) = (j) + 27ri?],^^'^^''** in the low-energy limit. 
Nowadays, the terminology "TLL" means the low-energy 
properties of such fermion models or the corresponding 
free boson theories. The c = 1 CFT consists of all the 
free boson theories with an arbitrary radius. The radius 
of the massless Dirac fermion (A2) is l/\/47r- 

We mention the scaling dimensions x ^ Al + Ar 
and conformal spins s = Al — Ar of primary fields in 
the c = 1 CFT, where Al(R) is the left (right) confor- 
mal weight. In our notation, the conformal weights of 
the vertex operator exp[i{aL(t>L +OiR(t>R)] is {Al,Ar) = 
(■g;^,-g^). Possible values of cil and cir can be deter- 
mined by the modular invariance. The currents Jl and 
Jr always have weights (1, 0) and (0, 1), respectively. In 
the Dirac fermion (Al), ^i, and ViJ have weights (|,0) 
and (0, ^), respectively. 



1 



(A4a) 



APPENDIX B: ISING MODEL AND c = 1/2 CFT 

We review some facts in terms of the 2D statistical 
(or ID transverse) Ising model and c = 1/2 CFT. The 
latter field theory emerges as the effective field theory in 
Sec. IV. 



1. Ising model 

Here, we siimmarizc relations between the Ising model 
on the square lattice and its continuum limit. (The con- 
tents here almost follows Ref. 116.) 

It is well known that both the order-disorder phase 
transition in the 2D Ising model and the quantum phase 
transition in the ID transverse Ising modcl*^^'*^*^ belong 
to the universality of the c = 1/2 CFT.38.39.ii9-i2i xhese 
two models are connected with each other via the transfer 
matrix method. The latter Hamiltonian is 



(Bl) 



where a" is a component of Pauli matrices settling in site 

J, and h is the transverse field. The critical point lies in 
h = J: \i {) < h < J {h > J), the order parameter sat- 
isfies ((7|) 7^ (=0). Here, let us introduce the disorder 
operator /Uj_|_i/2 on the dual lattice {j -1-1/2} as 



(B2a) 



p=i 

i-i 

'^i = n ^1 = Mj-i/2M^-M/2> (B2b) 

where A*J_fi/2 obey the same commutation relations as 
Pauli matrices. From Eqs. (B2), we have 

nTi[{cT]};J,h] = WTi[K+i/2};/i,J]. (B3) 

This is called the Kramers- Wannier duality. The self- 
dual point is just the critical one h = J . The ordered 
phase in the original Ising model, where ((t|) ^ 0, corre- 
sponds to the disordered phase in the dual model, where 
{^^^jJ^l|^2) — 0, and vice-versa. In the vicinity of the criti- 
cal point, continuum limits of C7| and A*^_)_i/2 respec- 
tively, associated with the order field <j{z, z) and disor- 
der field /i(z, z) in the c = 1/2 CFT, which are identical 
with Ua{z,z) and ^a{z,z) respectively in Sec. IV. It is 
worth while emphasizing that the relation between c7| 
and Mfe+i/2 nonlocal: these two commute when j > k, 
but anticommute otherwise. 

Real fermion operators can be introduced as 



Tlj = f^iMj-1/2. 



(B4) 
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where {'nj,Vk} = {CjXk} = 25jk, {vjXk} = 0. Inversely 
Ising operators arc written by fcrmions as 

= IJ'j+i/2 = -KjVj+i, (B5a) 

The Hamiltonian (Bl) can be described by these 
fermions, and is solvable in the fermion language. In 
particular, considering the vicinity of the critical point 
h = J, one can transform it to the following field- 
theoretical Hamiltonian: 



H 



TI 



where v = 2Jao, m = 2{h — J), x = j X ao and oq is the 
lattice constant. New Majorana (real) fermions and 
are defined as 



a = + 



2, = {C-V 



(B7) 



where r]{x) = rij/y/2ao and ({x) = (j/y/2ao. At the crit- 
ical point, the mass term vanishes, and the eflcective field 
theory (B6) then becomes a massless Majorana fermion 
model, which is just a c = 1/2 CFT. The fields £,l and 
are corresponding to ^£ and in Sec. IV B, respectively. 
From Eqs. (B5b) and (B7), it is obvious that ^l{r) has a 
nonlocal relation with a and /z, too. 

Two copies of the critical Majorana fermion theories 
are equivalent to the massless Dirac fermion (Al), i.e., 
a c = 1/2 + 1/2 = 1 CFT. Here let us denote the fields 
of two c = 1/2 systems, the corresponding Dirac fermion 
and boson theories as {^i^R,(7i, Hi), {£.l,r^(^2, H2), i'L,R 
defined as Eq. (37), and ^l,r, respectively. The U(l) 
currents in Eq. (A4) are written by Majorana fermions 
as follows: 



Energy operators £1,2 = i^L^£.R^ are mapped to 

£1 + £2 = i{lpLlpR - V'JjV'i) 

= I cos(v47r0), 

7ra 



(B8) 



(B9) 



where r\L,R are Klein factors in Eq. (A3). In addition, it 
is believed that order and disorder fields are bosonized 

g^g38,116,122 

cri(T2 ~ sin(\/7r(/)), ^1^2 ^ cos(V7r(/)), 

(Ti^2 cos(-\/7r0), /iicr2 ~ sin(V7''0), (BIO) 

where d is the dual of 6. 



the left (right) mover ^h(K) of the Majorana fermion, the 
energy operator e(z, z) = i^L£,R, the order field a and the 
disorder one /i. These are derivcd'^^'^^^ by making use of 
fusion rules and the Abelian bosonization based on the 
fact that two copies of c = 1/2 CFTs form a massless 
Dirac fermion model, a c = 1 CFT. The results are 

a(^)a(0) ~ x^, (Blla) 



27rz' 



(Bllb) 



aiz,z)aiO,0) ~ + ^\zf/^e{0,0), (BUc) 



H{z,z)n{0,0) 
a{z,z)fx{0,0) 

H{z,z)a{0,0) 



7r|z|3/4e(0,0), (Blld) 



|^|l/4 

^[eW4,V2^40) 



1^1 V4' 

+e-'^/^z'^^^R{0)], (Bile) 



t/4^1/2 



a(o) 



1^1 V4' 

+e-/*^V2^^(0)], (Bllf) 
(Bllg) 



a(^)a(o,o) ~ 


gi7r/4 


m(o,o) 


20Fzi/2 


U{z)aiO,0) ~ 




Mo,o) 


20F2^/2 


a(2)M(0,0) ~ 


g-i7r/4 


a(0,0) 




^Riz)i,{0,0) ~ 


gi7r/4 


a(0,0) 


20Fzi/2 



(Bllh) 
(Blli) 
(Bllj) 



As a reflection of nonlocal natures among Cl{r), and 
/Lt, some OPEs have a branch cut. The OPE (Bile) in- 
dicates that the product of the order and disorder fields 
must have a fermionic property. Following Ref. 112, in 
the main text, we often use a rule that a disorder field 
anticommutes with other disorder fields (when we con- 
sider some copies of c = 1/2 CFTs) and fermion fields, 
but commutes with itself and order fields. This is respon- 
sible for the improvement (33). 

As mentioned in Sec. IV, three copies of c = 1/2 CFTs 
are equivalent to the level-2 SU(2) WZNW model. Using 
the definition of the SU(2) currents (31), OPEs (B2a) and 
(B2b), we obtain the following OPEs among the currents: 



Jl{z)Jl{^) 

JI{z)j'r{q) 



Sab 



Sab 
477^^2 



+ ieabc 



2ttz ' 
2nz 



(B12) 
(B13) 



OPE 



APPENDIX C: SU(2) SYMMETRY IN SPIN 
SYSTEMS AND THE WZNW MODEL 



We write down the OPEs among the primary fields 
in the c = 1/2 CFT:i^'39'ii2,i22 ^^le identity operator 1, 



We consider what transformation for the effective field 
theories (28) or (35) is corresponding to the global spin 
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rotational transformation for spin-1 systems. 

The latter can be represented by a vector rotation form 

as 

Sj ^ T,i-ip3)Ty{-ip2)TA-^-i)S,=nSj, (CI) 

where Sj = '^{S^,S^,S^) and Ta{ip) stands for the 3D 
rotation about a axis by angle ifi [an SO (3) matrix]. For 
example, is defined as 



of the WZNW model. Reversely, both terms in Eq. (34) 
violate the chiral symmetry, and are invariant just under 
the "diagonal" rotation. 

Next, we focus on the rotation of the spin staggered 
part (C3b). When one represents the action of the 
WZNW model by using the matrix field gr^38,39 global 
SU(2) symmetry means that the action is invariant under 



(C5) 



cos if sin ip 
Tz{ip) = I — sin<^ costj? 
1 



(C2) 



From Eq. (32), the spin rotation on the lattice is inter- 
preted, in the WZNW model, as 



Jl,r 
I 



(C3a) 
(C3b) 



where J„ = ^(J^JlJ^) and / = ^(^1,52,^3) « 
'^(o'iA*2M3; o'3MiM2)- The transformation for the 
spin uniform part (C3a) is reproduced by the following 
transformation for Majorana fermions: 



L,R, 



(C4) 



where = "^(CiiCai^a)- One can confirm that the ef- 
fective theory for the Heiscnbcrg chain [Eq. (28) plus 
Eq. (34)] is invariant under the rotation (C4). Espe- 
cially, the level-2 SU(2) WZNW model (28) is invariant 
both under the rotation of left movers ^ and 
under that of right movers T^^R- These two sym- 

metries must correspond to the chiral SU(2) symmetry 



where F is a SU(2) matrix. A natural expectation is 

that the transformation (C5) leads to the rotation (C3b). 
Following this idea, in fact, one can verify that if the 
matrix V is parametrized as follows: 



V = 



-le 



»(v3-yi)/2 



sm 



y2 



-le 

gi(vP3+Vi)/2 cos^' 



(C6) 



the explicit correspondence between (C3b) and (C5) ap- 
pears. From a rotation (C3b), one also finds the spin 
rotation (CI) does not affect Qi = (J\<J2<J3- Proving that 
the action (28) and operators Oj in Table III are invari- 
ant under the spin rotation is an easy work. 

We touch the rotation around spin z axis Sj — > 
Tz{~(p)Sj. It does not affect the z component of the spin, 
and provides the transformations ^;!^ — > cos tp£,]^ — sin tp^^ 
and sin ipQ^+cos ipQ^. In the Dirac fcrmion picture, 
defined by Eq. (38), these U(l) transformations are 



i}L,R exp{iip)ipL,R- 



(C7) 



This is used in Sec. IV C. 



APPENDIX D: EFFECTIVE HAMILTONIAN IN SEC. IV B 



The effective Hamiltonian for the spin-1 ladder (1) with a uniform field (2a) under the condition H > m is 



Weff = [HireeiL, R, ^a, V') + ^free(i, R, H ^) + HLt{L, R, L, R) + HI,{L, R, L, R; ) 

+HUL, R, L, R; i;, ^) + HU^^, ; V>, i,) + H^i^, ^) + 4tag(a, n) 



(Dl) 



where 



Hiree = iv' (L^ O^L - R^ d^R) + ^j^- —d^ + m 



2m 



H)i> + zx(?ia.?i - eRd.£,l) + ^mCie|, 



(D2a) 
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Hl^^ = -^UlU^{{R^R-RR^f + {L^L-LL^f + {R^Rf + {L^Lf} 

-^{Ul + Ut){{R^R - RR^){L^L - LL^) + (i? ^ R){L L)} 

-XUlul{R^RLL^ + L^LRR^ + {R,L^ R, L)} 

^^l^{lJl + Ulf{{R^R - RR^){R -yR) + {L^L - W){L ^ L)} 

+^^{Ul + Ulf{{R^R - RR^){Dl - W) + {R^ R){L L)} 

+AJ_LaoUlUl [RL^W + LR^RL% (D2b) 
Hl^ = iXU+U.{[R^R-RR^+L^L-LL^]eLeR + [R^R + L^mli%} 

+J±ao{[Ul{R^R + RR^) + Ul{L^L + LL^MrH + [Ul{R^R + RR^) + Ul{L^L + LDMUI} 
-iJ±aoU+U-{[R^R - RR^ + L^L - lD]£,U\ + [R)R - RR^ + Dl - LL^iUn}, (D2c) 

HL = ^{Ul + Ul){{i,^-^i)^)[R^R-RR^ + L^L-W] + {'ij;^'^)[R^R + L^L]} 

o 

+^U+U_{il}'<^{R'^R + L^L) + ^p^p\RR^ + LL'<) + {R, L,^; ^ R, L, xjj)} 

-il^(lJl + f/!){(V'V - ^pil^^)[R^R - RR^ + L^L - W] + {R, L,ip^R, L, ij)} 

+^^{U+ - U-fii^^^R'^V + it^t) + ^^[R^R + L^L) + h.c.}, (D2d) 

- ^'mui + eji) - i^i^'i' + i'^miii + iUD, (D2e) 

hL = -^{(VV-V'W + (V'-V;)} + ^(V'V-V'V'^)(^-V;), (D2f) 

■ffstag = C^J±ao fi;«;ao/ia+i/Xa+2a'oMa+l/ia+2- (D2g) 



We define a new velocity v' = {H"^ — m^yl'^vlB.. In 
the interaction terms ii\~<^ , Dirac fermions L and R 
stand for L{—x) and R{—x), respectively. From the 
field-theoretical point of view, biquadratic terms such as 
{L^L — LLi^Y should be interpreted as a product of two 
point-splitting terms. For example, {L^ L — LL^)"^ means 
lims^o{L^L - LL^){-x) x {L^ L - LL^){-x + 5). 



APPENDIX E: ABELIAN BOSONIZATION IN 
THE SPIN-1 LADDER 



The low-cncrgy action for the spin-1 AF ladder is given 
in Eq. (35). If we make the boson theory with a scalar 
field (j)a from two Ising systems cja and da {a = 1, 2, 3), we 
can obtain a bosonized Hamiltonian for the action (35) 
through several relations in Appendices A and B. It has 



been already given in Ref. 17. The result is 

3 . 

Wlad = V / dx{'^[Ill + {d,<Paf]-—COs{V^cPa) 

a—l 

+Bi [ COs{V47T4>a) COs{V47T4>a+l) 

+ cos{V^0a) cos(y47r6'a+i)] 

+B2[lla'n-a+l + dx4>adx<f>a+l] 

-I--B3 [ sin ( V4^(j)a ) sin ( \/47r0a+ 1 ) 

+ sm{VI^0a) sm{VI^0a+i)] (El) 

+B4, sin(i/7r(/)a) cos(y7r0a+i) C0s(v^(?!)a+2)}, 

where Si oc j^^, ^2,3 oc ^ and B4 ~ CfJjLao. 
We can not consider the Klein factors correctly, because 
the formulas (BIO) are not perfect. However, here let us 
assume that it is allowed. Actually, the same Hamilto- 
nian in Ref 17 (though its notation is diff'erent from ours) 
surely seems to work well. 

Like Ref. 17, we take the following mean-field prescrip- 
tion for Wiad- (i) We neglect Bi,2,3 terms, and leave only 
the free part (fermion kinetic and mass terms) and the 
most relevant interaction B^ term, (ii) For the weak AF 
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rung-coupling case, the mass potential — cos(\/47r(^a) will 
be still dominant, and therefore the field (pa is locked near 
the point = 0. This argument allows the following ap- 
proximation: cos(Y^0a) B {& constant). Through (i) 
and (ii), the Hamiltonian Tiiad is reduced to 



where 11$ and Hi, are canonical conjugated momenta of 
and respectively. The most relevant rung-coupling 
term J±aQC^KK(TalJ.a+ifJ'a+2d-afia+ifi'a+2 is mapped as 
follows: 



njMF 



dx{^-[Ill + id,<l>,f] (E2) 

a— 1 

-— C0s(\/4^^a) + B^BiSm{y/^(j)a)]. 

This is three copies of a double-sine-Gordon model. 

For details of this mean- field argument, see Ref. 17. 

Next, we consider the case that three boson fields 
4>a, 5" and $ are composed of {(Ja-,(Ja), (ca+i, o'a+2) and 
((Ta+i, (7a+2), respectively. Like Eq. (El), the fermion 
free part in the total Hamiltonian is mapped to 

Wfree = j dx^ [(H^ + (d^C^af) + (£1% + {d,^f) 



+(n?+(a.$)^) 



m 



cos(v47r^a) 



cos(v47r<i?) -I- cos(v47r<l?) 



(E3) 



J±[sm{^/^T4>a) cos(\/7r$) cos(\/7r$) 
-|- cos{^/^^^a){ cos{^/nQ) cos{^/nQ) 
+ sin(x/7re) sin(\/7re)}] 
JlB cos [x/7r(e-e)]. 



(E4) 



On the right side of the arrow — we performed 

the same mean-field approximation used in Eq. (E2): 
cos{y/TT(pa) B and sm.{^/^T(t)a) 0. 
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